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Introduction  and  identification  of  the  challenges 

There  has  been  a  recent  explosion  of  interest  in  using  terahertz  (THz)  radiation  to 
various  applications  in  chemistry,  biology,  physics  (spectroscopy),  medicine,  material 
science  and  imaging.  This  interest  is  motivated  by  the  fact  that  the  THz  part  of  the 
spectrum  is  energetically  equivalent  to  many  physical,  chemical  and  biological  processes. 
Simultaneously  there  has  been  a  rapid  development  of  a  broad  array  of  experimental  tools 
for  working  with  THz  radiation.  With  the  advent  of  modem  micro-fabrication  and 
MEMS  technologies,  the  operating  frequencies  of  vacuum  electron  microwave  devices 
are  being  extended  into  the  THz  regions.  One  example,  a  1  THz  Backward  Wave 
Oscillator  (BWO)  is  commercially  available1.  However,  state-of-the-art  THz  BWOs  are 
inefficient  (0.001%)  and  very  heavy  (>350kg,  including  required  magnets  and  power 
supplies),  making  them  impractical  for  widespread  applications. 

Extending  microwave  sources  into  the  sub-mm  wave  region,  combined  with  the 
need  for  miniaturization,  poses  a  variety  of  very  serious  challenges  for  electron  cathode 
performance,  as  well  as  for  electron  beam  transport,  quality  and  focusing.  Extending  the 
frequency  to  the  THz  range  requires  confining  the  electron  beam  to  an  increasingly  small 
cross  section.  In  fact,  since  the  required  beam  current  density  is  proportional  to  the 
operating  frequency  squared  ,  the  challenge  is  to  generate  and  realize  in  the  interaction 
space  the  extremely  high  current  densities  (~10  A/cm  )  required  even  for  low  power 
(-0.1W)  sources.  A  second  challenge  is  to  overcome  the  traditional  reliance  on  bulky 
magnets  and  solenoids  needed  for  magnetic  confinement  of  electron  beams. 

Achieving  the  desired  electron  beam  current  density  needed  for  compact  vacuum 
THz  sources  is  very  challenging.  Explosive  emission  cathodes  offer  high  current  density, 
but  suffer  from  short  pulse  length  (typically  <lpsec),  poor  beam  quality  and  fast 
deterioration  of  electrodes.  Conventional  thermionic  cathodes  widely  used  in  the 
microwave  tube  industry  can  provide  long  pulse  operation,  but  typically  are  limited  to 
low  current  density  (<  10  A/cm2). 
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In  this  research  program  we  investigated  a  novel  approach  to  fill  the  need  for 
advanced,  miniaturized,  plasma-assisted  electron  sources  suitable  for  tomorrow’s  THz 
microwave  vacuum  (MVE)  systems.  To  establish  the  context,  let  us  note  that  there  are 
some  alternative  concepts  of  electron  guns  proposed  for  THz  sources.  The  approach 
investigated  in  this  research  program  includes  two  significant  innovations: 

1)  the  use  of  a  miniaturized,  long  lived,  high  current  density,  plasma- 
assisted  electron  sources,  and 

2)  the  accompanying  use  of  ion  focusing  to  eliminate  the  need  for  heavy 
focusing  magnets  and  improve  the  quality  of  the  electron  beam. 

Plasma-assisted  cathodes  look  attractive  because  of  their  ability  to  produce  high 
current  densities.  They  can  be  considered  “equivalent”  to  indestructible  thermionic 
cathodes  operating  at  effective  temperature  ~20,000°C  having  a  zero  work  function 


Plasma  density  [cm'^]  X  Temperature172  [eV] 

Fig.  1.  Plasma-assisted  cathodes  are  characterized  by  very 
high  current  density  capability,  as  needed  for  THz  sources. 

without  perveance  limits,  that  do  not  exist  in  the  plasma.  For  example,  Fig.  1  illustrates 
the  calculated  maximum  current  density  from  plasma-assisted  cathodes  under  various 
conditions.  This  figure  describes  a  hypothetical  situation  where  all  available  electrons 
are  extracted  from  the  plasma.  A  limit  of  200  A/cm2  corresponds  to  a  1 .5- 1013  cm'3,  2eV 
plasma. 
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Plasma-assisted  cathodes  can  be  compatible  with  high-vacuum  technology 
necessary  for  THz  sources;  a  short  discussion  of  this  topic  is  in  order.  The  “common 
wisdom”  is  that  plasma-assisted  cathodes  are  fundamentally  incompatible  with  extreme 
vacuum  conditions  because  they  require  a  gas  pressure  ~10'4  Torr,  which  is  prohibitively 
high  for  conventional  vacuum  devices.  Note,  however,  that  practically  all  vacuum 
microwave  devices  contain  some  amount  of  plasma.  So  the  issue  is  therefore  not  the 
plasma  itself,  but  rather  how  much  plasma  is  acceptable.  As  long  as  the  plasma 
frequency  of  the  plasma  in  the  interaction  space  of  a  device  is  much  less  than  the 
operating  frequency  (cop  « co  or,  more  exactly,  co2«co2),  the  device  can  be 

practically  considered  as  a  vacuum  device.  Since  in  the  THz  frequency  range  the 
operating  frequency  is  much  higher  than  the  frequency  of  conventional  L-band  or  S-band 
microwave  devices,  this  condition  can  be  easily  satisfied  in  THz  frequency  band,  even  in 
the  presence  of  a  background  gas  with  the  pressure  on  the  order  of  10'4  Torr,  as  shown  in 
Table  1. 

Table  1:  THz  MVE  devices  driven  by  plasma-assisted  cathodes  can  be  compatible  with 

vacuum  electronic  device  technology. 


GHz  MVE  devices 

driven  by  thermionic 

cathodes 

THz  MVE  devices 

driven  by  plasma- 

assisted  cathodes 

Background  gas  pressure 

~10'8  Torr 

~10‘4  Torr 

Background  gas  density 

310s  cm'3 

Background  plasma  frequency 

(assuming  full  ionization) 

~108Hz 

~101U  Hz 

Ratio  of  plasma  frequency  to 

operating  frequency 

<10'2 
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As  a  result,  THz  MVE  sources  are  much  less  sensitive  to  the  presence  of 
background  plasma  than  GHz  MVE  sources.  For  this  reason,  THz  MVE  devices  driven 
by  plasma-assisted  electron  cathodes  are  projected  to  be  fully  compatible  with  practical 
vacuum  device  concepts  mandated  by  the  MVE  community.  In  addition,  plasma-assisted 
cathodes  are  compatible  with  high  pulse  repetition  rates,  and  even  with  continuous  wave 
(CW)  operation.  They  are  scalable,  can  be  miniaturized,  require  neither  a  filament,  nor  a 
filament  power-supply,  and  therefore  can  be  turned  on  instantly  without  any  warm-up 
time.  They  can  possibly  be  gated  using  simple,  low-voltage  grids,  compatible  with 
moderate  vacuum  technology  and  are  relatively  insensitive  to  impurities.  Finally,  no 
externally  applied  magnetic  field  is  needed,  as  an  ion  channel  produced  by  beam  impact 
ionization  of  the  background  gas  can  be  used  to  guide  and  transport  the  beam  through  the 
interaction  region.  For  these  reasons,  plasma-assisted  cathodes  seem  to  be  attractive 
candidates  to  drive  vacuum  THz  sources. 

Finally,  miniaturized  plasma-assisted  cathodes  could  prove  to  be  a  breakthrough 
enabling  technology  for  THz  MVE  devices,  with  the  potential  to  advance  the 
technological  and  scientific  base  and  to  impact  military,  commercial,  industrial  and 
scientific  applications  through  the  development  of  transferable,  commercially  viable 
technologies.  For  example,  present  THz  sources  are  currently  notoriously  inefficient 
(1/100  of  1%).  With  our  new  approach,  there  is  an  exciting  opportunity  for  efficiency 
enhancement  of  THz  sources  by  at  least  two  orders  of  magnitude  (or  more)  because  of 
the  extremely  high  current  density  and  high  perveance  capabilities  of  plasma-assisted 
cathodes. 

Most  significant  accomplishments 

The  most  significant  accomplishments  are  summarized  below.  (For  more  details,  see 
Appendix  A,  and  for  selected  publications  see  Appendix  B.) 

Experimental  accomplishments 

The  most  significant  experimental  accomplishment  was  the  demonstration  that 
plasma-assisted  miniature  cathodes  can  generate  high  perveance  electron  beams  at 
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relatively  low  voltages  (>1000A/cm2,  0.5-3kV),  and  that  those  beams  could  propagate 
without  external  guiding  magnetic  fields  over  distances  compatible  with  the  length  of  the 
RF  structure  of  THz  sources.  This  enables  the  design  of  future  THz  sources  operating 
with  relatively  high  efficiency  at  high  power  levels. 

1)  Two  types  of  miniature,  plasma-assisted  cathodes  were  designed,  fabricated  and  tested 
during  this  research  period.  Some  of  the  most  important  characteristics  of  these  guns  are 
summarized  in  Table  2.  Note  that  very  high  current  densities  were  successfully 
demonstrated,  -1000  [A/cm2],  as  shown  in  Table  2  and  in  Fig.  2. 

Table  2:  The  most  important  characteristics  of  plasma-assisted  electron  guns  studied 

during  research  period 


Beam  voltage  [kV] 

0.5-3 

Beam  current  [A] 

0-0.3 

Maximum  current  density  [A/cm2] 

-1000 

2)  Beam  propagation  over  a  distance  of  10-20  mm  was  demonstrated  in  the  absence  of 
any  external  guiding  magnetic  field.  These  distances  are  compatible  with  the  interaction 
length  of  THz  devices,  so  this  demonstration  is  very  encouraging. 

One  type  of  a  plasma-assisted  cathode  was  described  in  a  conference  article  3. 

For  more  information  on  both  versions  of  plasma  cathodes  see  Appendix  A. 
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Fig. 2.  Demonstrated  current  density  as  the  function  of  time.  Extracted 
electron  current  (red)  and  the  plasma  discharge  current  (black)  in  the  glass 
capillary  are  almost  equal.  (For  details  see  Appendix  A.) 
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Theoretical  accomplishments 

In  the  framework  of  the  MiPRI  program  our  group  had  experimentally 
demonstrated  some  possibilities  to  use  miniature  plasma  electron  guns  for  generating 
electron  beams  of  high  current  density.  In  view  of  these  accomplishments,  the  theoretical 
efforts  were  focused  on  analyzing  the  possibilities  to  generate  coherent  radiation  in  the 
THz  frequency  range  with  a  reasonably  high  efficiency  and  at  high  power  levels.  Our 
analysis  was  restricted  by  compact  devices  capable  of  operation  at  low  voltages.  Since  at 
frequencies  approaching  1  THz  it  is  much  simpler  (if  not  the  only  possible)  to  fabricate 
planar  microwave  circuits,  we  considered  devices  having  a  planar  geometry,  which  most 
conveniently  can  be  driven  by  sheet  electron  beams.  Since  high-power  operation  at  very 
short  wavelengths  is  possible  without  breakdown  only  in  the  case  when  the  interaction 
space  is  characterized  by  dimensions  much  larger  than  the  wavelength,  our  analysis  was 
focused  on  devices  using  highly  selective,  oversized  resonators  (like  the  orotron  whose 
schematic  is  shown  in  Fig.  3). 


Output 


Electron  Beam 


Fig.  3.  Schematic  of  the  orotron:  an  open  resonator  formed  by  a  concave  mirror  and  a 
periodic  grating  plate  is  excited  by  a  sheet  electron  beam. 


1.  The  interaction  region  in  the  devices  specified  above  has  transversely  nonuniform 
interaction  space.  This  nonuniformity  should  decrease  the  interaction  efficiency. 
To  evaluate  the  importance  of  this  effect  the  analytical  study  had  been  performed. 
The  analysis  performed  allowed  us  to  formulate  requirements  to  design 
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parameters,  under  which  the  effect  of  this  nonuniformity  on  device  efficiency  is 
relatively  weak.  This  work  was  reported 4  and  published  5. 

2.  After  a  certain  analysis  it  became  clear  that  the  device  performance  is  severely 
limited  by  the  fact  that  the  field  of  a  slow  wave  moving  with  the  phase  velocity 
close  to  electron  velocity  is  localized  near  the  slow-wave  circuit  in  a  layer  whose 
thickness  is  much  smaller  than  the  wavelength.  To  position  a  sheet  electron  beam 
within  this  layer  and  avoid  the  beam  interception  with  the  circuit  requires  very 
fine  beam  alignment  and  strong  guiding  magnetic  field.  It  also  requires  very  high 
current  density  for  providing  in  such  a  narrow  electron  beam  enough  current  for 
exciting  microwave  oscillations  in  a  circuit.  To  avoid  this  limitation,  two 
possibilities  had  been  studied. 

3.  First,  we  considered  the  case  when  the  beam  current  exceeds  the  start  current  only 
slightly  that  results  in  excitation  of  microwave  oscillations  with  small  amplitude. 
In  such  a  case,  interaction  with  the  slow  wave  causes  small  changes  in  electron 
energy.  However,  this  small  modulation  of  electron  energies  by  the  RF  field  can 
result,  at  long  enough  distances,  in  significant  electron  bunching  improving  the 
device  efficiency.  Corresponding  analytical  theory  allowing  us  to  make  some 
estimates  for  optimal  operating  parameters  of  the  orotron-type  devices  was 
developed;  results  were  published  6  and  will  be  presented  as  an  invited  keynote  at 
the  upcoming  31st  Int.  Conf.  on  IR&MM  Waves  &14th  Int.  Conf.  on  THz  Electronics, 
Shanghai,  P.R.  China  1 .  It  was  found,  for  example,  that  in  the  orotron  operating  at 
the  frequency  of  375  GHz  and  driven  by  a  4kV,  0.3  A  electron  beam  having 
electron  current  density  of  300  A/cm2  it  is  possible  to  produce  up  to  13  W  output 
radiated  power  with  the  efficiency  exceeding  1%.  Note  that  such  a  density  can  be 
realized  either  in  non-convergent  electron  beams  generated  by  cathodes  with  high 
current  density  (such  as  described  above  plasma  cathodes  developed  at  the 
University  of  Maryland)  or  in  convergent  electron  beams  generated  by  more 
conventional  cathodes  with  relatively  low  current  density.  Also  note  that  a 
relatively  low  value  of  the  output  efficiency  is  due  to  the  fact  that  significant  part 
of  the  microwave  power  radiated  by  an  electron  beam  should  be  lost  due  to  ohmic 
losses  in  the  structure  walls  (the  situation  typical  for  so  high  frequencies). 
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4.  Secondly,  we  analyzed  the  case  when  an  electron  beam  is  slightly  inclined  with 
respect  to  the  grating  plate.  The  device  of  such  a  configuration  is  known  as  the 
clinotron  (or  “klinotron”).  In  such  a  configuration  shown  in  Fig.  4,  when  the 
device  length  is  long  enough,  all  electrons  of  an  arbitrarily  thick  electron  beam 
sooner  or  later  enter  the  layer  occupied  by  a  slow  wave  near  the  grating  plate. 
This  fact  allows  one  to  greatly  increase  the  beam  cross-section,  mitigate 
requirements  to  the  electron  current  density  and  beam  alignment  and  realize  much 
higher  power  levels  of  operation.  (Of  course,  in  the  case  of  continuous  wave 
operation  some  limitations  may  be  imposed  by  the  density  of  the  beam  power 
deposited  to  the  grating  plate  in  the  area  of  beam  interception.)  The  theory 
describing  operation  of  such  devices  has  been  developed  recently.  The  paper 
describing  it  is  in  the  process  of  preparation  for  publishing  in  Physics  of  Plasmas. 


Fig.4.  Schematic  of  the  clinotron:  a  thick  electron  beam  propagates  at  a  small  angle  with 
respect  to  the  grating  plate  and  interacts  with  a  slow  wave  localized  near  the  grating 
before  hitting  the  plate. 
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Appendix  A 

Two  types  of  plasma-assisted  miniature  electron  guns  for  THz 

vacuum  devices. 
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Abstract 

Efforts  to  improve  efficiency  in  terahertz  vacuum  electronic 
sources  such  as  backward-wave  oscillators  (BWO)  have 
motivated  studies  of  cathodes  capable  of  generating  high- 
density  electron  beams  so  that  various  interaction  circuits  may 
be  driven  closer  to  saturation.  Preliminary  results  from 
experiments  on  a  miniature  plasma  cathode  demonstrate  that 
current  densities  of  up  to  1  kA/cm2  in  a  100-micron  diameter 
beam  have  been  achieved.  Initial  analysis  of  a  850-GHz, 
slotted  slow-wave  structure  show  that  the  circuit  could  be 
driven  at  3-4  times  the  starting  current  and  that  efficiencies  of  a 
few  percent  could  be  realized.  Tests  on  a  gridded  plasma 
cathode  show  that  stable  beams  of  200  mA  at  a  pulse  length  of 
100  psec  and  repetition  frequencies  of  up  to  200  Hz  were 
generated. 

Introduction 

Terahertz  sources,  such  as  the  clinotron  and  orotron,  are  based 
on  Cerenkov  emission  from  the  interaction  between  an  electron 
beam  and  the  electromagnetic  field  in  a  slow  wave  structure 
(SWS)  within  an  open  resonator  [1].  However,  typical 
efficiencies  even  in  the  most  advanced  devices  are  in  the  104  - 
10'3  range.  This  is  due,  in  part,  because  saturated  operation  is 
difficult  to  achieve  since  devices  with  thermionic  cathodes  are 
limited  in  terms  of  current  density.  Large-area  plasma  cathodes 
which  were  developed  for  high-power  helix  BWO’s  have  been 
demonstrated  to  reliably  generate  high  current  densities  (50- 
500  A/cm2)  [2],  In  the  present  work,  increased  output  power 
and  improved  efficiency  in  terahertz  sources  will  be  pursued  by 
developing  miniaturized  plasma  cathodes  capable  of  driving 
oscillators  at  several  times  the  starting  current.  To  illustrate,  a 
simple  calculation  of  the  starting  current  in  an  open-resonator 
SWS  follows. 


Calculation  of  Starting  Current 


In  [3],  it  was  found  that  the  ratio  of  the  amplitude  of  the 
synchronous  space  harmonic  of  the  wave  to  the  amplitude  of 
the  zero  space  harmonic  is  maximal,  when  the  depth  of  SWS 
grooves  is  equal  to  a  quarter  of  wavelength  and  the  width  of 
grooves  is  equal  to  the  half  of  the  structure  period.  The 
minimum  starting  current  in  the  device  with  a  gaussian  axial 
structure  of  the  operating  mode  can  be  determined  according  to 
[3]  as 


I  — *h-cRy 

Ast  ^  Vb  ’ 

xCohm  int 


(i) 


Here  Sb  is  the  beam  cross-sectional  area,  L,„,  is  the  interaction 
length,  p0  is  the  velocity  of  the  electrons  normalized  to  the 
speed  of  light,  c,  Vh  is  the  beam  voltage  and  Qohm  is  the  ohmic 
quality  factor.  It  has  been  assumed  that  total  Q  is  dominated  by 
ohmic  losses,  which  yields  the  following  formula: 


Qohm  8skO  +  ?Jd)’  (2) 

where  H  is  the  distance  between  the  two  mirrors  forming  the 
resonator,  Ssk  is  the  skin  depth,  which  for  copper  at  850  GHz 
is  about  0.15  microns,  X  is  the  radiation  wavelength  and  d  is  the 
period  of  the  SWS.  Note  that  the  ratio  of  the  wavelength  to  the 
structure  period  can  be  expressed,  consistent  with  the  condition 
of  Cherenkov  synchronism,  in  terms  of  the  electron  velocity. 
With  H  =  2mm,  the  ohmic  Q-factor,  depending  on  the  choice 
of  voltage  and  grating  period,  is  in  the  range  from  1 .3  *  103  to 
1.7-  103.  Thus  the  minimum  start  current  density  can  be 
determined  as 

(l  +  3/30)Vb(kV)-  0) 

H  Ant(£"2M) 

It  follows  from  (3)  that  if  the  separation  between  mirrors  is  3 
mm  and  the  interaction  length  is  5  mm,  the  start  current  density 
for  a  850  GHz  clinotron  operating  10  kV  is  equal  to  64  A/cm  . 
For  high  efficiency  operation  the  current  density  should  exceed 
the  start  current  density  by  at  least  a  factor  of  2,  which  results 
in  the  range  of  nominal  current  densities  of  about  130  A/cm2. 
For  the  case  of  the  beam  width  of  3  mm  and  beam  thickness  of 
0. 1  mm,  these  numbers  correspond  to  the  total  beam  current  of 
200-300  mA. 

Simple  Plasma  Cathode 

The  plasma  cathode,  which  is  shown  schematically  in 
Fig.  1,  was  a  scaled  version  of  the  design  given  for  high- 
power  applications  referenced  in  [2].  The  hollow  cathode 
assembly  was  fabricated  from  oxygen-free  copper  with  a 
beam  aperture  of  1  mm  in  diameter  and  a  gas  port  in  the 
back  wall  in  which  a  2  mm  capillary  tube  was  soldered. 


Fig.  1 :  Expanded  cross-sectional  diagram  of  cathode  pieces 


The  grid  was  fabricated  from  thin  molybdenum  sheet 
and  was  electrically  isolated  from  the  cathode  by  a  mica 
insulating  disk.  Both  had  100  micron  holes  for  beam 
extraction  that  were  aligned  concentrically  with  the  axis 
of  the  hollow  cathode  during  assembly.  The  assembled 
gun  is  shown  in  Fig.  2  mounted  in  a  vacuum  test 
chamber. 


Fig.  2:  Miniature  plasma  cathode  shown  mounted  in  a 
vacuum  gun  test  chamber 


The  grid  was  driven  by  a  solid  state  pulse  generator  at  voltages 
in  the  range  of  500-800  V,  which  was  determined  to  be 
adequate  to  ionize  the  fill  gas  in  the  hollow  cathode  in  less  than 
0.5  psec  and  pulse  lengths  in  the  range  10-200  psec.  The 
current  in  the  grid-cathode  discharge  was  found  to  be 
approximately  twice  the  beam  current.  The  gun  was  tested 
according  to  the  parameters  given  in  Table  I. 

Table  I.  Experimental  Parameters  for  Gun  Test 


Beam  Voltage 
Beam  Current 


0.5- 1.1  kV 
0  -  0.3  A 


Max  Current  Density 

Plasma  Current 
Gas  (constant  source) 
Cathode  Pressure 
Base  Pressure 
Pumping 


~  1  kA  /  cm^ 
0.5  A 
N,  Ar,  Xe 
10-50  microns 
<  0.01  microns 
Turbo  w/diff. 


A  schematic  of  the  test  circuit  is  shown  in  Fig.  3  along  with 
typical  voltage  and  current  waveforms  in  Fig.  4.  The  voltage 
was  measured  between  the  cathode  and  anode  giving  the  total 
accelerating  potential  including  the  DC  bias  applied  between 
cathode  and  anode.  Figure  5  is  a  plot  of  the  beam  current 
measured  by  a  wideband  current  monitor  installed  around  a 
ground  return  lead  connected  to  the  anode  structure.  The 
diameter  of  the  beam  was  measured  to  be  approximately  70 
microns  by  a  witness  plate  located  2  mm  downstream  of  the 
extraction  aperture. 


remained  stable  for  many  months  of  operation  over  a  wide 
range  of  operating  conditions. 

Hollow  Cathode  Grid 


Fig.  3:  Schematic  of  gun  test  circuit  which  shows  the 
approximate  pressures  in  the  hollow  cathode  and 
beam  transport  regions 
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Fig.  5:  Measured  current-voltage  characteristics  of  plasma 
gun 
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At  this  writing,  the  gun  has  been  operated  for  periods  of  up  to 
12  hours  with  pulse  lengths  up  to  200  psec  and  repetition 
frequencies  up  to  100  Hz,  corresponding  to  a  maximum  duty 
factor  of  2%  with  only  radiation  cooling.  During  the  first  few 
minutes  of  operation,  the  ionization  characteristics  of  newly 
fabricated  guns  drifted  somewhat,  and  the  gas  pressure  had  to 
be  decreased  to  maintain  the  desired  current  and  voltage.  After 
this  short  “conditioning”  period,  the  gun  characteristics 
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Abstract 

Efforts  to  improve  efficiency  in  terahertz  vacuum  electronic 
sources  such  as  backward-wave  oscillators  (BWO)  have 
motivated  studies  of  cathodes  capable  of  generating  high- 
density  electron  beams  so  that  various  interaction  circuits  may 
be  driven  closer  to  saturation.  Preliminary  results  from 
experiments  on  a  miniature  plasma  cathode  demonstrate  that 
current  densities  of  up  to  1  kA/cm2  in  a  100-micron  diameter 
beam  have  been  achieved.  Initial  analysis  of  a  850-GHz, 
slotted  slow-wave  structure  show  that  the  circuit  could  be 
driven  at  3-4  times  the  starting  current  and  that  efficiencies  of  a 
few  percent  could  be  realized.  Tests  on  a  gridded  plasma 
cathode  show  that  stable  beams  of  200  mA  at  a  pulse  length  of 
100  psec  and  repetition  frequencies  of  up  to  200  Hz  were 
generated. 

Introduction 

Terahertz  sources,  such  as  the  clinotron  and  orotron,  are  based 
on  Cerenkov  emission  from  the  interaction  between  an  electron 
beam  and  the  electromagnetic  field  in  a  slow  wave  structure 
(SWS)  within  an  open  resonator  [1],  However,  typical 
efficiencies  even  in  the  most  advanced  devices  are  in  the  10"4  - 
10'3  range.  This  is  due,  in  part,  because  saturated  operation  is 
difficult  to  achieve  since  devices  with  thermionic  cathodes  are 
limited  in  terms  of  current  density.  Large-area  plasma  cathodes 
which  were  developed  for  high-power  helix  BWO’s  have  been 
demonstrated  to  reliably  generate  high  current  densities  (50- 
500  A/cm2)  [2].  In  the  present  work,  increased  output  power 
and  improved  efficiency  in  terahertz  sources  will  be  pursued  by 
developing  miniaturized  plasma  cathodes  capable  of  driving 
oscillators  at  several  times  the  starting  current.  To  illustrate,  a 
simple  calculation  of  the  starting  current  in  an  open-resonator 
SWS  follows. 


Calculation  of  Starting  Current 


In  [3],  it  was  found  that  the  ratio  of  the  amplitude  of  the 
synchronous  space  harmonic  of  the  wave  to  the  amplitude  of 
the  zero  space  harmonic  is  maximal,  when  the  depth  of  SWS 
grooves  is  equal  to  a  quarter  of  wavelength  and  the  width  of 
grooves  is  equal  to  the  half  of  the  structure  period.  The 
minimum  starting  current  in  the  device  with  a  gaussian  axial 
structure  of  the  operating  mode  can  be  determined  according  to 
[3]  as 


I  =— — ^h-cRV 
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Here  5),  is  the  beam  cross-sectional  area,  Lin,  is  the  interaction 
length,  Po  is  the  velocity  of  the  electrons  normalized  to  the 
speed  of  light,  c,  Vb  is  the  beam  voltage  and  Qohm  is  the  ohmic 
quality  factor.  It  has  been  assumed  that  total  Q  is  dominated  by 
ohmic  losses,  which  yields  the  following  formula: 


Qohm  Ssk(3  +  A  /  d)  ’  (2) 

where  H  is  the  distance  between  the  two  mirrors  forming  the 
resonator,  Ssk  is  the  skin  depth,  which  for  copper  at  850  GHz 
is  about  0.15  microns,  X  is  the  radiation  wavelength  and  d  is  the 
period  of  the  SWS.  Note  that  the  ratio  of  the  wavelength  to  the 
structure  period  can  be  expressed,  consistent  with  the  condition 
of  Cherenkov  synchronism,  in  terms  of  the  electron  velocity. 
With  H  =  2mm,  the  ohmic  Q-factor,  depending  on  the  choice 
of  voltage  and  grating  period,  is  in  the  range  from  1.3  •  103  to 
1.7  •  103.  Thus  the  minimum  start  current  density  can  be 
determined  as 


Js,  (A/cm2)  =  2-1  o4  (1  +  3  /?„ )  Vb(kV)  ■  (3) 

H  Lmt(cm ) 

It  follows  from  (3)  that  if  the  separation  between  mirrors  is  3 
mm  and  the  interaction  length  is  5  mm,  the  start  current  density 
for  a  850  GHz  clinotron  operating  10  kV  is  equal  to  64  A/cm2 
For  high  efficiency  operation  the  current  density  should  exceed 
the  start  current  density  by  at  least  a  factor  of  2,  which  results 
in  the  range  of  nominal  current  densities  of  about  130  A/cm2. 
For  the  case  of  the  beam  width  of  3  mm  and  beam  thickness  of 
0.1  mm,  these  numbers  correspond  to  the  total  beam  current  of 
200-300  mA. 


Simple  Plasma  Cathode 

The  plasma  cathode,  which  is  shown  schematically  in 
Fig.  1,  was  a  scaled  version  of  the  design  given  for  high- 
power  applications  referenced  in  [2].  The  hollow  cathode 
assembly  was  fabricated  from  oxygen-free  copper  with  a 
beam  aperture  of  1  mm  in  diameter  and  a  gas  port  in  the 
back  wall  in  which  a  2  mm  capillary  tube  was  soldered. 


Fig.  1 :  Expanded  cross-sectional  diagram  of  cathode  pieces 


The  grid  was  fabricated  from  thin  molybdenum  sheet 
and  was  electrically  isolated  from  the  cathode  by  a  mica 
insulating  disk.  Both  had  100  micron  holes  for  beam 
extraction  that  were  aligned  concentrically  with  the  axis 
of  the  hollow  cathode  during  assembly.  The  assembled 
gun  is  shown  in  Fig.  2  mounted  in  a  vacuum  test 
chamber. 


Fig.  2:  Miniature  plasma  cathode  shown  mounted  in  a 
vacuum  gun  test  chamber 


The  grid  was  driven  by  a  solid  state  pulse  generator  at  voltages 
in  the  range  of  500-800  V,  which  was  determined  to  be 
adequate  to  ionize  the  fill  gas  in  the  hollow  cathode  in  less  than 
0.5  psec  and  pulse  lengths  in  the  range  10-200  psec.  The 
current  in  the  grid-cathode  discharge  was  found  to  be 
approximately  twice  the  beam  current.  The  gun  was  tested 
according  to  the  parameters  given  in  Table  I. 


Table  I.  Experimental  Parameters  for  Gun  Test 


Beam  Voltage 
Beam  Current 


0.5- 1.1  kV 
0-0.3  A 


Max  Current  Density 

Plasma  Current 
Gas  (constant  source) 
Cathode  Pressure 
Base  Pressure 
Pumping 


~  1  kA / cm^ 
0.5  A 
N,  At,  Xe 
10  -  50  microns 
<  0.01  microns 
Turbo  w/  diff. 


A  schematic  of  the  test  circuit  is  shown  in  Fig.  3  along  with 
typical  voltage  and  current  waveforms  in  Fig.  4.  The  voltage 
was  measured  between  the  cathode  and  anode  giving  the  total 
accelerating  potential  including  the  DC  bias  applied  between 
cathode  and  anode.  Figure  5  is  a  plot  of  the  beam  current 
measured  by  a  wideband  current  monitor  installed  around  a 
ground  return  lead  connected  to  the  anode  structure.  The 
diameter  of  the  beam  was  measured  to  be  approximately  70 
microns  by  a  witness  plate  located  2  mm  downstream  of  the 
extraction  aperture. 


remained  stable  for  many  months  of  operation  over  a  wide 
range  of  operating  conditions. 
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Fig.  3:  Schematic  of  gun  test  circuit  which  shows  the 
approximate  pressures  in  the  hollow  cathode  and 
beam  transport  regions 
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Fig.  5:  Measured  current-voltage  characteristics  of  plasma 
gun 
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At  this  writing,  the  gun  has  been  operated  for  periods  of  up  to 
12  hours  with  pulse  lengths  up  to  200  psec  and  repetition 
frequencies  up  to  100  Hz,  corresponding  to  a  maximum  duty 
factor  of  2%  with  only  radiation  cooling.  During  the  first  few 
minutes  of  operation,  the  ionization  characteristics  of  newly 
fabricated  guns  drifted  somewhat,  and  the  gas  pressure  had  to 
be  decreased  to  maintain  the  desired  current  and  voltage.  After 
this  short  “conditioning”  period,  the  gun  characteristics 
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To  analyze  such  microwave  sources,  2.5D  and  3D  codes  are 
usually  used. 
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-Analysis  of  the  effect  of  transverse  nonuniformity  of  the  RF 
field  on  the  efficiency  of  microwave  sources  (this  presentation) 


Formulation  of  the  problem 

A  beam  of  electrons  moving  linearly  along  the  device  axis 
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The  function  Ezs  describes  transverse  nonuniformity  of  the 
RF  field  acting  on  electrons 


Formulation  of  the  problem  (cont.) 
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What  can  we  do  with  this  general  formula? 


Quasi-linear  approximation 
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Quasi-linear  theory  (cont.) 
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nonuniformity  of  the  RF  field  on  the  interaction  efficiency 


EXAMPLES 


RESULTS 


Maximum  value  is  71 


RESULTS  (cont.) 


RESULTS  (cont.) 
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The  normalized  thickness  2.0  corresponds  to  the  real  beam  thickness  of  0.022mm. 


RESULTS  (cont.) 
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RESULTS  (cont.) 
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This  paper  contains  a  simple  analytical  theory  that  allows  one  to  evaluate  the  effect  of  transverse 
nonuniformity  of  the  rf  field  on  the  interaction  efficiency  in  various  microwave  sources  driven  by 
linear  electron  beams.  The  theory  is,  first,  applied  to  the  systems  where  the  beams  of  cylindrical 
symmetry  interact  with  rf  fields  of  microwave  circuits  having  Cartesian  geometry.  Also,  various 
kinds  of  microwave  devices  driven  by  sheet  electron  beams  (orotrons,  clinotrons)  are  considered. 
The  theory  can  be  used  for  evaluating  the  efficiency  of  novel  sources  of  coherent  terahertz 
radiation.  ©  2005  American  Institute  of  Physics.  [DOI:  10.1063/1.2034427] 


I.  INTRODUCTION 

Many  important  sources  of  coherent  microwave  radia¬ 
tion  (such  as  klystrons,  traveling-wave  tubes,  backward- 
wave  oscillators,  magnetrons,  and  gyrotrons),  at  least  in  their 
simplest  configurations,  have  a  cylindrical  symmetry  of  the 
interaction  space.  This  fact  greatly  simplifies  the  theory  of 
such  devices  and,  in  many  cases,  allows  developers  to  do  at 
the  first  step  of  design  one-dimensional  (ID)  simulations. 
Then,  some  additional  factors,  such  as  transverse  motion  of 
electrons  moving  initially  linear  along  the  device  axis,  asym¬ 
metry  caused  by  the  presence  of  input  and  output  couplers, 
etc.,  can  be  taken  into  account. 

In  many  cases,  however,  new  demands  for  escalating 
performance  characteristics  of  the  existing  devices  as  well  as 
an  interest  to  devices  based  on  new  principles  of  operation 
lead  to  analysis  of  configurations,  in  which  the  interaction 
space  does  not  have  an  axial  symmetry.  Orotrons1,2  as  well 
as  sheet-beam3,4  klystrons  can  serve  as  examples  of  such 
configurations.  As  an  example  of  demands  resulting  in  the 
preference  of  devices  without  axial  symmetry  of  the  interac¬ 
tion  space,  one  can  mention  an  existing  strong  interest  in  the 
development  of  sources  of  coherent  terahertz  radiation.  At  so 
high  frequencies  it  is  much  simpler  (if  not  the  only  doable)  to 
fabricate  planar  microwave  circuits  (see,  e.g.,  Ref.  5)  instead 
of  cylindrical  ones. 

In  such  configurations  of  microwave  sources,  where 
electron  beams  moving  initially  along  the  device  axis  inter¬ 
act  with  the  electromagnetic  fields  having  different  ampli¬ 
tude  at  location  of  different  beamlets,  the  effect  of  this  trans¬ 
verse  nonuniformity  on  the  device  efficiency  is  usually 
calculated  with  the  use  of  two-and-one-half-dimensional6 
(2.5D)  and  three-dimensional7  (3D)  codes.  For  initial  evalu¬ 
ation  of  the  effect  of  this  nonuniformity  on  the  device  effi¬ 
ciency,  however,  some  analytical  methods  can  be  proposed. 
These  methods  should  be  treated  as  complementary  to  3D 
code  simulations.  The  purpose  of  such  methods  is  to  help 
designers  to  easily  evaluate  the  effect  of  transverse  nonuni¬ 
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formity  of  the  rf  field  on  the  device  efficiency  before  starting 
some  time-consuming  3D  simulations. 

In  this  paper,  we  propose  some  of  such  methods  and  give 
some  estimates  of  the  effect  of  transverse  nonuniformity  on 
the  efficiency  in  several  device  configurations.  The  paper  is 
organized  as  follows.  In  Sec.  II,  the  formalism  describing  the 
quasilinear  theory  is  presented.  Such  a  theory  was  originally 
used  for  evaluating  the  effect  of  transverse  nonuniformity  of 
the  rf  field  on  the  efficiency  of  the  gyrotron.8  In  Sec.  Ill,  a 
number  of  examples  are  considered  that  illustrate  the  appli¬ 
cation  of  the  method  to  some  specific  configurations.  In  Sec. 
IV,  a  simple  nonlinear  theory  describing  similar  effects  in  the 
clinotron  is  developed.  In  Sec.  V  we  apply  the  developed 
formalism  to  the  design  of  a  W-band  traveling-wave  tube 
(TWT),  which  was  done  with  the  use  of  ID  and  2.5D 
simulations,10  and  compare  the  results  of  our  analytical  treat¬ 
ment  with  the  data  of  numerical  simulations.  Finally,  Sec.  VI 
contains  a  brief  discussion  and  summary. 


II.  QUASILINEAR  THEORY 
A.  General  formalism 

Let  us  consider  a  beam  of  electrons  moving  linearly 
along  the  device  axis  with  the  velocity  vz  and  interacting 
with  a  synchronous  space  harmonic  of  the  axial  electric  field 
of  a  resonator, 

Ez  =  ReiAEjrJMe^'-W}.  (1) 

Here  the  functions  f(z )  and  Esz( rx)  describe  the  axial  and 
transverse  structures  of  this  field,  respectively. 

Electron  motion  in  such  a  ID  case  can  be  described  by 
equations  for  the  electron  energy  e  and  time  that  a  particle 
spends  in  the  interaction  space: 
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Below  we  will  consider  a  stationary  electron  beam,  for 
which  in  (2)  d/dt=d/ dt+vz(dl dz)  =  vz(d/ dz).  Hence,  intro¬ 
ducing  a  slowly  variable  phase  difference  between  electrons 
and  the  phase  of  the  synchronous  wave,  and  nor¬ 

malized  variables  z'  =  (co/c)z,  y  =e/mc2,  f3=vz/c,  h=kzc/o), 
and  A'=eA/mc<a,  we  can  rewrite  (2)  and  (3)  omitting  primes 
as 


II 

^  1 

-R s{AEJeiS}, 

(4) 

df) 

y 

(5) 

dz 

~r=  h. 

V72-  1 

In  (5),  the  general  relation  between  electron  energy  and  ve¬ 
locity,  y=  1  /  \j\-01,  was  used.  Equations  (4)  and  (5)  should 
be  supplemented  by  the  initial  conditions:  y(0)  =  y0  and 
r)(0)  =  r/0>  where  the  entrance  phase  ri0  is  uniformly  distrib¬ 
uted  from  0  to  2 v.  Below  we  will  assume  that  the  beam  is 
monoenergetic,  i.e.,  the  spread  in  the  initial  energies  is  neg¬ 
ligibly  small. 

Transverse  nonuniformity  of  the  rf  field,  which  was  dis¬ 
cussed  in  the  Introduction,  is  described  in  Eq.  (4)  by  the 
function  Esz(r±).  Correspondingly,  the  efficiency  of  a  beam, 
whose  density  distribution  in  the  device  cross  section  can  be 
given  by  the  function  'l'(r^),  can  be  determined  as 


Note  that  the  function  ‘'f/(rx)  obeys  the  normalization  condi¬ 
tion  Ssb'^r(r±)dsx=Sb,  where  Sb  is  the  cross-sectional  area  of 
the  beam.  Also  in  (6)  L=wL/c  is  the  normalized  length  of 
the  interaction  space. 


B.  Quasilinear  approximation 


terms  at  right-hand  sides  proportional  to  A  are  much  smaller 
than  the  first  ones.  These  perturbations  in  electron  motion,  as 
follows  from  (4)  and  (5),  are  equal  to 

y(1)  =  R  e{eifl°A£»,  (8) 


#(I)  (/0-l)3/2i 


Re{e,l\4Ev,v}. 


(9) 


Here  we  introduced  u=p{)felSz  dz.’  and  v=fz(lu(z')dz'.  First- 
order  perturbations  in  electron  energy,  being  substituted  in 
the  formula  for  the  device  efficiency  (6),  yield  zero  after 
averaging  over  the  entrance  phases.  Therefore,  even  for  ob¬ 
taining  the  region  of  parameters  where  the  self-excitation  of 
the  device  is  possible,  it  is  necessary  to  calculate  the  second- 
order  terms,  i.e.,  y{2)  proportional  to  A2,  as  it  was  done  in 
Ref.  11.  Representing  the  electron  energy  as  7=7(o)-7(i) 
-  y(2)  and  linearizing  the  right-hand  side  of  Eq.  (4)  up  to  the 
second-order  terms,  one  can  readily  obtain  the  solution  for 


To  calculate  the  efficiency,  one  should  average  this  change  in 
energy  over  the  entrance  phases  and  transverse  coordinates 
of  beamlets.  Resulting  expression  for  the  efficiency  valid  in 
the  framework  of  the  linear  theory  is  the  following: 


’7lin  =  |A|2" 


2(y0-i)U-i) 


3/2 


CRe  J-\ 


lin  ■ 


(10) 


Here 


C=f  [  ^(rx)|£J,(r1)|2^1  (11) 

bbJsh 

is  the  coefficient  describing  the  coupling  of  an  electron  beam 
to  the  synchronous  component  of  the  rf  field  and  Re  Tym  is 
the  real  part  of 

Jrun  =  -i[  f(z)e~iSzv(z)dz,  (12) 

Jo 


In  the  framework  of  a  quasilinear  theory,  we  will  con¬ 
sider  the  effect  of  the  rf  field  on  electron  motion  by  a  method 
of  successive  approximations  assuming  the  amplitude  of  the 
rf  field,  A,  in  Eqs.  (4)  and  (5)  to  be  a  small  parameter.  Note 
that  this  assumption,  as  follows  from  normalization  of  the 
amplitude,  corresponds  to  the  restriction  on  the  original  am¬ 
plitude  of  the  rf  field  A(kV/cm)<t2TT-51 1/X(cm),  which  is 
practically  always  fulfilled.  The  first  steps  in  this  derivation 
are  quite  similar  to  those  described  in  Ref.  11. 

In  the  zero-order  approximation,  electron  energy  and 
phase  are  equal  to 

7(0)  =  7o>  #(0)  =  #o+  &■  O) 

Here  we  introduced  initial  mismatch  of  synchronism  be¬ 
tween  electrons  and  the  wave,  S=l/f30-h. 

First-order  terms  in  electron  motion  follow  from  repre¬ 
sentations  7=  7(o)- 7(i)  and  #=#(o)+#(1),  in  which  the  last 


which  is  the  complex  “linear”  gain  function  equivalent  to  the 
beam-loading  conductance  known  in  the  theory  of  klystrons 
(see,  e.g.,  Ref.  12). 

Next,  we  shall  calculate  the  saturation  terms  in  the  de¬ 
pendence  of  the  efficiency  on  the  rf  field  amplitude.  For  do¬ 
ing  this,  it  is  necessary  to  move  with  the  method  of  succes¬ 
sive  iterations  further  up  to  the  fourth-order  term  in  electron 
energy,  because  the  third-order  terms  yield  zero  after  averag¬ 
ing  over  the  entrance  phase,  as  it  happened  with  the  first- 
order  terms  considered  above.  Omitting  the  details  of  rather 
lengthy  derivation,  we  will  give  here  the  final  expression  for 
the  saturation  terms  in  the  efficiency,  which  in  the  framework 
of  the  quasilinear  theory  can  be  presented  as 

7=%n-%at-  (13) 

For  a  device  with  long  enough  interaction  length  ( L=  coLIc 
>  1),  this  saturation  term  ?/sat  is  equal  to 
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%at  =  H4 


1 


8(ro-D(^-D9/2 
Here  the  coefficient 


D  Re 


(14) 


d  =  t[  ^(ri)l£«(ri)|4^i  (15) 

6 bJsb 

describes  the  effect  of  transverse  nonuniformity  of  the  rf 
field  on  the  saturation,  and  Re  ,Fsat  is  the  real  part  of  the  gain 
function  describing  the  saturation  in  the  framework  of  the 
quasilinear  theory: 


Note  that  this  formula  (16)  coincides  with  the  definition  of 
the  self-saturation  coefficient  derived  in  the  theory  of  mode 
interaction  in  gyrotrons.13 

Substituting  all  these  expressions  into  (13)  results  in  the 
following  formula  for  the  efficiency  in  the  quasilinear  ap¬ 
proximation: 


1 


’7=2(y0-l)(^-l)34HC^in 

1 


4(^-1) 


M\4DT'S  J. 


(17) 


Here  the  primes  denote  the  real  parts  of  the  gain  functions, 
and  the  coefficients  C  and  D  determined  above  characterize 
the  effect  of  the  transverse  nonuniformity  of  the  rf  field  on 
the  device  efficiency  in  the  small-signal  regime  and  in  the 
saturation  stage,  respectively.  From  (17)  one  can  easily  find 
the  optimal  intensity  of  the  rf  field,  M=\A\2: 

WoPt  =  2(^-l)3^^.  (18) 

•'"sat 


Maximum  efficiency  that  corresponds  to  this  intensity  is 
equal  to 


U-i  )3/2TO2c 
2(ro-D  Kt  o' 


At  low  voltages,  the  first  term  at  the  right-hand  side  of  (19) 
can  be  reduced  to  (-)^-l)3/2/2(y0-l)  —  [2(eVywc2)]I/2.  As 
follows  from  (19),  the  effect  of  the  transverse  nonuniformity 
of  the  rf  field  on  the  device  maximum  efficiency  can  be 
estimated  by  the  ratio 


(20) 


So,  the  maximum  efficiency  depends  on  the  difference  be¬ 
tween  the  effects  of  transverse  nonuniformity  on  the  linear 
and  nonlinear  properties  of  an  oscillator,  which  are  described 


TABLE  I.  Typical  configurations  of  microwave  sources  with  transversely 
nonuniform  interaction  space. 


by  the  coefficients  C  and  D,  respectively.  The  application  of 
this  formula  to  various  specific  configurations  of  microwave 
circuits  and  electron-beam  cross  sections  is  given  in  Sec.  III. 

111.  EXAMPLES 

In  this  and  next  sections  we  will  study  the  effect  of  the 
transverse  nonuniformity  of  the  rf  field  in  the  device  configu¬ 
rations  that  are  schematically  shown  in  Table  I.  The  first  one 
is  a  cylindrical  pencil-like  electron  beam  inserted  in  the  rf 
structure  formed  by  two  parallel  plates.  (Such  a  structure  can 
model  the  cavities  used  in  sheet-beam  klystrons.)  The  second 
example  is  the  same  beam  inserted  in  a  rectangular  cavity. 
The  third  example  is  a  sheet  electron  beam  propagating  over 
the  surface  of  grating  plate.  The  case  of  parallel  propagation 
of  this  beam,  which  models  the  interaction  in  the  orotron,1,2 
will  be  considered  in  this  section.  The  case  of  a  tilted  beam, 
which  models  the  clinotron9’14  and  needs  a  certain  modifica¬ 
tion  of  the  theory,  will  be  considered  in  Sec.  IV. 

A.  Cylindrical  beam  between  two  parallel  plates 

As  well  as  in  open  Fabry-Perot  resonators  used  in  lasers, 
the  rf  field  in  the  wide  transverse  (y)  dimension  can  be  de¬ 
scribed  by  Hermite  polynomials.  Below  we  will  restrict  our 
consideration  by  the  lowest-order  Hermite  polynomial, 
which  is  the  Gaussian  profile,  because  the  modes  with  a 
larger  number  of  variations  in  y  direction  have  higher  dif¬ 
fraction  losses  through  the  open  ends.  The  field  distribution 
between  the  plates,  i.e.,  in  x  direction,  can  be  represented  as 
a  standing-wave  pattern.  Below  we  will  consider  a  mode 
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*  with  one  variation  in  this  direction,  because  modes  with  a 
larger  number  of  such  variations  have  much  higher  frequen¬ 
cies.  Thus,  the  dependence  of  the  axial  component  of  the 
electric  field  will  be  given  as 

£,  =  e"^")2  cosf“j,  (21) 

where  the  parameters  a  and  b  characterize  the  scale  of  the  rf 
field  nonuniformity  in  y  and  x  directions,  respectively,  and  jc 
varies  from  -b  to  b  between  the  plates.  Regarding  an  elec¬ 
tron  beam  we  assume  that  this  is  a  cylindrical  solid  beam  of 
a  radius  R. 

Substituting  (21)  and  the  beam  coordinates  into  (11)  and 
(15)  yields  the  following  expressions  for  the  coefficients  C 
and  D,  respectively: 

C=  —  f  e~2p“'^\  \jl  -i?  +  —  sin(p6y'  l  -  i?)  \dv,  (22) 
Wo  [  Ph  J 

D  =  ■ —  f  e~4p“pl\  3  \ll  -  i/2  +  —  sinfp/jVl  -  i/2) 

2W0  {  Pb 

+  —  sin(2pil\/l  -  i?)  \dv.  (23) 

2  Pb  J 

In  (22)  and  (23),  we  introduced  normalized  variables  and 
parameters:  v=y/R,  pa=R/a,  and  ph=irR/b.  Then,  expres¬ 
sions  given  by  (22)  and  (23)  were  substituted  into  Eq.  (20). 

The  results  of  the  calculations  are  shown  in  Fig.  1.  Here 
Fig.  1(a)  shows  the  dependence  of  the  factor  / j,  which  char¬ 
acterizes  the  effect  of  the  transverse  nonuniformity  on  the 
device  efficiency  for  our  model,  on  the  parameter  pb  for  sev¬ 
eral  values  of  pa,  and  Fig.  1(b)  shows  the  contours  of  equal 
values  of  I\  in  the  plane  of  parameters  pb  vs  pa.  Note  that  the 
maximum  value  pi  max=rr  corresponds  to  the  beam  diameter 
equal  to  the  distance  between  the  plates.  (Let  us  call  the 
readers’  attention  to  the  difference  in  notations  for  pa  and  ph 
by  a  factor  of  i r.)  As  follows  from  these  data,  as  far  as  the 
beam  radius  does  not  exceed  2a,  which  is  the  width  of  the 
Gaussian  distribution  at  the  e~]  level,  and  also  is  not  larger 
than  one-third  of  the  distance  between  the  plates,  the  effi¬ 
ciency  decreases  by  less  than  a  factor  of  2. 


B.  Cylindrical  beam  in  a  rectangular  cavity 


The  axial  component  of  the  rf  electric  field  in  a  rectan¬ 
gular  cavity  for  the  case  of  odd  number  of  field  variations  in 
transverse  dimensions  can  be  given  as 


(24) 


(In  the  case  of  even  m  or  n  cos  should  be  replaced  by  sin.) 
Note  that  Eq.  (24)  is  written  for  the  reference  frame  with  the 
origin  in  the  center  of  a  cavity. 

Assuming  that  the  beam  center  coincides  with  this  ori¬ 
gin,  one  can  substitute  (24)  into  (11)  and  (15)  and  using  a 
given  geometry  of  a  pencil-like  beam  readily  arrive  at  the 
following  expressions  for  coefficients  C  and  D: 


Beam  radius  normalized  to  the  horizontal  nonuniformity,  PB 

FIG.  1.  Effect  of  nonuniformity  on  the  interaction  efficiency  of  a  cylindrical 
electron  beam  located  between  two  parallel  plates:  (a)  the  factor  Jj  as  the 
function  of  the  parameter  pb  characterizing  the  field  nonuniformity  between 
two  plates  for  several  values  of  the  parameter  pa  characterizing  the  field 
variation  in  a  wide  transverse  direction;  (b)  contours  of  constant  values  of/, 
in  the  plane  of  parameters  pb  vs  p„. 


1  1  ,  ,  1  ,  , 

C-  +  •iri(Po)+  Jliph) 
4  2  pfl  2  pb 


+  —  |  cos(p„v)sin(piVl  -  i?)dv,  (25) 

^PbJ  0 

D= - f  [cos(2pai)  +  4  cos(pflv)  +  3]  X  I  3  Vl  -  v2 

16ttJ0  [ 

4  _  i  _  I 

+  —  sin(piyl  -vL)  +  —  sin(2pfcVl  -  W  f dv.  (26) 

Pb  2Pb  J 

Here  pa=mrR/a  and  ph=nnrR/b.  In  (25)  we  used  the  known 
integral  (see,  e.g.,  3.752.2  in  Ref.  15):  -u2 cos(pv)dv 

=  (Tr/2p)7,(p). 

Substituting  these  formulas  for  C  and  D  into  Eq.  (20) 
readily  yields  the  value  of  the  ratio  I2  characterizing  the  ef¬ 
fect  of  transverse  nonuniformity  on  the  efficiency.  Contours 
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FIG.  2.  Effect  of  nonuniformity  of  the  interaction  efficiency  of  a  cylindrical 
electron  beam  in  a  rectangular  cavity:  contours  of  constant  values  of  the 
factor  /2  are  shown  in  the  plane  of  parameters  characterizing  the  field  non¬ 
uniformity  in  both  transverse  directions. 


of  /2= const  are  shown  in  Fig.  2  in  the  plane  pb  vs  p„.  Of 
course,  these  contours  are  symmetric  with  respect  to  pb=pa- 
Recall  that  for  modes  with  one  variation  the  maximum  value 
of  a  corresponding  parameter  p  does  not  exceed  tt.  Within 
pa  h  *£  7 t  the  factor  I2  is  larger  than  0.7,  i.e.,  efficiency  degra¬ 
dation  is  not  very  strong.  Note  that  the  dependence  of  I2  on 
the  parameter  ph  is  essentially  the  same  as  the  dependence  of 
/j  shown  in  Fig.  1  because  this  normalized  variable  pb  has 
the  same  meaning  in  both  cases.  At  the  same  time,  the  ratio 
I2  is  less  sensitive  to  the  variable  pa  than  the  ratio  /,  because 
of  the  difference  in  normalization  of  the  beam  radius  in  this 
variable  to  the  horizontal  scale  of  the  rf  field  nonuniformity 
in  both  cases.  For  instance,  if  we  assume  that  the  scale  of 
horizontal  nonuniformity,  a'A\  of  the  Gaussian  distribution  in 
case  A  is  equal  to  the  distance  between  sidewalls,  aiB\  in 
case  B  and  the  rf  field  in  the  latter  case  has  only  one  variation 
in  y  direction,  we  get  the  following  relation  between  vari¬ 
ables  pa  for  both  cases:  p(B>=Trp(Al.  So,  in  this  specific  ex¬ 
ample  the  values  of  I2  shown  in  Fig.  2  for,  let  say,  pa=7r 
should  be  compared  with  the  values  of  /,  shown  in  Fig.  1(b) 
for  pa=  1.  Of  course,  this  coefficient  of  proportionality  de¬ 
pends  on  the  position  of  sidewalls  in  case  B. 

C.  The  orotron 

In  microwave  sources  based  on  the  use  of  coherent  Cher¬ 
enkov  radiation,  electrons  interact  with  a  synchronous  elec¬ 
tromagnetic  wave  (or  synchronous  space  harmonic  of  it  in  a 
periodic  slow-wave  structure),  which  propagates  in  the  direc¬ 
tion  of  electron  flow  with  the  phase  velocity,  vph=w/kz, 
smaller  than  the  speed  of  light  c.  Correspondingly,  the  axial 
wave  number  kz  of  such  of  a  wave  is  larger  than  ai/c  and, 
hence,  such  a  wave  has  an  imaginary  transverse  wave  num¬ 
ber,  which  characterizes  the  decay  of  this  wave  with  the  dis¬ 
tance  from  the  surface  of  a  slow-wave  structure.  As  follows 
from  the  condition  of  Cherenkov  synchronism  between  this 
slow  wave  and  electrons,  the  size  of  the  region  of  transverse 
localization  of  such  a  slow  wave  can  be  estimated  as  (see, 


e.g„  Ref.  16)  L±  ~(y^-l)1/2X.  At  low  voltages,  where 
-\yn ~  (2eVbl me1)'12 <  1,  and  at  short  wavelengths,  this 
distance  is  very  small.  Therefore  the  effect  of  the  finite  thick¬ 
ness  of  an  electron  beam  becomes  well  pronounced. 

To  estimate  the  role  of  this  effect  on  the  interaction  ef¬ 
ficiency,  consider  a  simple  model,  in  which  only  the  effect  of 
the  transverse  exponential  decay  of  the  rf  field  of  a  slow 
wave  is  taken  into  account.  Such  a  model  is  shown  in  Table 
I.  First,  consider  the  case  of  the  orotron  where  an  electron 
beam  propagates  parallel  to  the  grating  plate. 

Let  us  represent  the  axial  rf  electric  field  as 

Ez  =  Ae~KX,  (27) 

where  k  is  the  absolute  value  of  the  wave  number  in  the  x 
direction.  Assume  that  a  sheet  electron  beam  has  a  finite 
thickness  A  and  its  midplane  is  located  at  the  distance  x* 
=*mjn+A/2,  where  xmin  is  the  clearance  between  the  beam 
and  the  grating  surface.  Edge  effects  in  the  y  direction  will 
be  neglected.  Consider  the  case  when  the  rf  field  amplitude  is 
optimal  for  electrons  located  in  the  midplane.  Then,  with  the 
use  of  (27)  one  can  readily  get  the  following  simple  formula 
for  the  factor  I: 


tanh(/rA) 
~  *A 


(28) 


As  follows  from  (28),  the  efficiency  degrades  by  a  factor  of 
2  when  the  normalized  thickness  kA  is  close  to  2.0.  For 
instance,  in  a  device  driven  by  a  5-kV  electron  beam  and 
operating  at  a  0.5-mm  wavelength,  this  normalized  thickness 
of  2.0  corresponds  to  the  real  beam  thickness  of  0.022  mm. 


IV.  THE  CLINOTRON 

As  is  known,  to  mitigate  the  issue  associated  with  effi¬ 
ciency  and  power  degradation  at  short  wavelengths  caused 
by  miniaturization  of  the  interaction  space  and  corresponding 
clearance  between  the  beam  and  the  grating  plate,  the  device 
called  the  clinotron  was  proposed.9  The  key  distinction  be¬ 
tween  the  clinotron  and  other  slow-wave  devices  (such  as 
backward-wave  oscillators  and  orotrons)  is  an  inclined  posi¬ 
tion  of  an  electron  beam  with  respect  to  the  surface  of  a 
slow-wave  structure.  In  such  a  device,  even  electrons  of  a 
relatively  thick  beam  can  sooner  or  later  interact  with  a 
strong  enough  rf  field  near  the  grating  surface.  As  shown  in 
Ref.  9  (see  also  Ref.  14  for  latest  references),  even  a  rather 
small  angle  of  inclination  (less  than  1  deg)  allows  one  to 
greatly  increase  the  radiated  microwave  power  and  facilitates 
the  device  alignment  and  the  requirements  to  the  focusing 
magnetic  field.  It  should  be  noted  that,  according  to  Ref.  9, 
first  papers  on  the  theory  and  experiments  with  clinotrons 
were  published  in  USSR  in  1958,  while  in  1957  Karp  pub¬ 
lished  his  paper  on  the  experiments  with  backward-wave  os¬ 
cillators  (BWOs),17  in  which  the  principle  of  the  clinotron 
and  corresponding  improvements  in  the  BWO  performance 
in  the  case  of  a  tilted  electron  beam  were  described. 

To  apply  our  formalism  to  the  clinotron,  introduce  the 
reference  frame  associated  with  a  beam  inclined  at  the  angle 
a  where  the  axial  z'  and  transverse  x'  coordinates  are  related 
to  the  coordinates  z  and  x  used  above  as: 
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z'  -z  cos  a  -x  sin  a,  x'  -z  sin  a  +  x  cos  a.  (29) 

At  small  angles  a,  the  most  important  effect  is  the  depen¬ 
dence  of  the  exponential  term  of  the  rf  field  acting  upon 
electrons  in  (27)  on  z':  E.=Ae~KX = Ae~~K(,:,l~az'\  Here  x0  is  the 
transverse  coordinate  of  an  electron  at  the  entrance. 

Let  us  now  modify  the  basic  Eqs.  (4)-(6)  for  the  case  of 
the  clinotron.  For  simplicity  assume  that  the  function  f(z)  in 
(4)  is  constant:  f(z)  - 1 .  Also  assume  that  perturbations  in 
electron  energy  are  small  [y=  y0-  7n),  where  1 7(1)1^  To], 
while  the  electron  bunching  caused  by  these  perturbations 
can  be  significant.  These  assumptions  reduce  Eqs.  (4)  and  (5) 
to 


<97(0  , 

-  Ae  x°+az  cos  if, 
dz! 

(30) 

<9fl  „  TO) 

*'  Vo-1  )3/2' 

(31) 

Here  we  normalized  x0  and  a  to  the  absolute  value  of  the 
transverse  wave  number  by  introducing  xQ=  kxq  and  a=xa 
[in  (30)  the  bars  are  omitted],  the  detuning  S  has  the  same 
meaning  as  above.  Correspondingly,  the  efficiency  (6)  can  be 
rewritten  as 

1  1  f'x0.roin+^  [  1  f2ir  j 

[ul  ™ 

Here  we  ignore  the  field  nonuniformity  in  the  wide  (y)  trans¬ 
verse  direction.  The  interaction  of  electrons  with  the  field 
ends  when  electrons  hit  the  grating  plate,  i.e.,  at  z^nd=x0/a. 
Thus,  electrons  with  different  transverse  entrance  coordi¬ 
nates  have  different  interaction  lengths. 

Equations  (30)  and  (31)  are  similar  to  known  equations 
describing  the  physics  of  interaction  between  electrons  and 
electromagnetic  field  in  such  devices  as  traveling-wave  tubes 
and  free-electron  lasers.  As  known,18’19  those  equations  for 
electron  motion  can  be  reduced  to  the  equation  for  the  non¬ 
linear  pendulum.  In  the  case  of  Eqs.  (30)  and  (31),  this 
means  introducing  a  new  set  of  variables:  y'  =  y(1j/ \A(y^ 
-1)3/4,  £=\(Az7(7o-l)3/4,  S'  =(y^-\)v‘iS/\^A,  and  a'  =  tf0 
- 1  )3/4a/  \A.  In  new  variables,  these  equations  do  not  contain 
two  parameters  (A  and  y0)  and  have  the  following  forms: 


2L  _  e-x0+a'€  ft 

(33) 

dd 

(34) 

S'  +  y’. 

<?£ 

Correspondingly,  the  efficiency  given  above  by  (32)  can  be 
rewritten  as 


V  = - —yAr/, 

7o-l 


(35) 


where  the  normalized  efficiency  rj  is  given  by 


These  simple  equations  are  derived  for  a  thick  electron  beam 
propagating  at  a  small  angle  with  respect  to  the  grating  plate, 
which  they  reach  at  fend=x0/a'. 

When  the  beam  thickness  is  negligibly  small  (all  elec¬ 
trons  have  the  same  transverse  entrance  coordinate  x0),  one 
can  reduce  these  equations  further  by  introducing  y=  y'e^1, 
g=£e~X:/2,  a=a'ex°n,  and  S=  S' ex°n.  This  change  of  vari¬ 
ables  reduces  (33)  and  (34)  to 


—  =  cos  ■&, 

(37) 

% 

dft  — 

—  =  S  +  y. 

(38) 

K 

Correspondingly,  the  efficiency  can  be  characterized  by  its 
normalized  value, 

1  f277  - 

V=~ ~J  (39) 

This  set  of  equations  contains  only  three  normalized  param¬ 
eters  that  are  proportional  to  the  interaction  length,  initial 
detuning  of  synchronism  between  electrons  and  a  slow  wave 
and  a  tilting  angle. 

The  results  of  the  calculations  for  the  clinotron  with  a 
thin  electron  beam  are  shown  in  Fig.  3  as  contours  of  equal 
values  of  the  normalized  efficiency  rj  in  the  plane  of  param¬ 
eters  “initial”  detuning  <5  versus  normalized  interaction 
length  £end  for  several  values  of  the  normalized  tilt  angle  a. 
Here  Fig.  3(a)  (a=0)  corresponds  to  the  orotron  with  a  beam 
parallel  to  the  grating  plate.  At  zero  detuning  or  mismatch  of 
synchronism  between  electrons  and  synchronous  wave,  the 
efficiency  is  zero  as  in  all  O-type  electron  devices.  This  ef¬ 
ficiency  is  positive  when  electron  velocity  slightly  exceeds 
the  phase  velocity  of  the  wave,  i.e.,  at  negative  d’s  [see  defi¬ 
nition  of  this  detuning  given  above  after  (7)]. 

Two  other  figures,  Figs.  3(b)  and  3(c),  characterize  the 
clinotron  operation  for  the  tilting  parameter  a  equal  to  0.1 
and  0.2,  respectively.  These  contours  clearly  demonstrate 
that  the  tilt  of  a  beam  results  in  significant  enhancement  of 
the  efficiency.  The  fact  that  this  normalized  efficiency,  which 
characterizes  the  efficiency  of  electron  bunching,  is  larger 
than  1 .0  is  not  important  because  the  efficiency  given  by  (35) 
contains  a  small  parameter  proportional  to  the  square  root  of 
the  rf  field  amplitude.  In  this  regard,  this  set  of  equations  is 
quite  similar  to  the  simplest  set  of  nonlinear  equations  for  the 
conventional  TWT  (Ref.  1 8)  where  the  efficiency  can  be  rep¬ 
resented  as  the  product  of  the  normalized  efficiency,  which 
can  be  larger  than  1 .0,  and  the  Pierce  gain  parameter,  which 
has  a  small  value;  see  also  pages  164  and  165  of  Ref.  20  with 
the  reference  to  Slater’s  “Microwave  Electronics.”  (Note 
that  the  relativistic  gyrotron  can  also  be  described  by  similar 
equations.21)  Coming  back  to  Fig.  3,  let  us  note  that  optimal 
parameters  and  the  maximum  normalized  efficiency  of  a  de- 
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FIG.  3.  Contours  of  constant  values  of  the  normalized  efficiency  in  the 
plane  of  parameters  “normalized  initial  detuning  vs  normalized  interaction 
length”  in  the  clinotron  with  different  values  of  the  tilt  parameter  5:  (a)  5 
=0,  (b)  5=0.1,  and  (c)  5=0.2. 


vice  with  a  parallel  beam  [Fig.  3(a),  5=0]  correspond  to  the 
results  shown  in  Fig.  2(a)  of  Ref.  22. 

Since  we  have  mentioned  above  that  the  rf  field  ampli¬ 
tude  plays  here  a  role  of  a  small  parameter,  it,  possibly, 
makes  sense  to  discuss  the  scaling  of  the  efficiency  and  ra- 


FIG.  4.  Electron  trajectories  in  the  phase  space  “electron  energy  vs  electron 
phase”  in  the  clinotron  with  5=0.1  (a),  and  5=0.2  (b);  the  normalized  initial 
detuning  and  the  normalized  interaction  length  have  optimal  values  in  ac¬ 
cordance  with  Fig.  3(c). 

diated  microwave  power  with  the  operating  voltage  in  more 
detail.  In  general,  the  balance  equation  for  the  power  of  mi¬ 
crowave  losses  and  the  power  withdrawn  by  the  rf  field  from 
an  electron  beam  can  be  written  as 

^W=  v  VbIb.  (40) 

Here  the  left-hand  side  represents  the  power  of  losses  of  the 
microwave  energy  W=j/l|2V’  stored  in  the  cavity  having  a 
volume  V  and  the  quality  factor  Q,  while  the  right-hand  side 
is  the  product  of  the  interaction  efficiency  and  the  beam 
power  VbIb.  When  the  efficiency  is  given  by  (35),  which  at 
low  voltages  reduces  to  t}rj:(A<l2IVlM)r/,  the  rf  field  ampli¬ 
tude,  in  accordance  with  (40),  is  proportional  to  A  *  V^2/^3. 
For  electron  guns  with  space-charge-limited  emission  this 
means  A  «  V^/2.  Correspondingly,  the  efficiency  (35)  scales  as 
r/xVV27/  and  the  radiated  microwave  power  scales  as  Frf 

*  Wb- 

In  Fig.  4,  the  electron  trajectories  are  shown  in  the  phase 
space  of  the  set  of  Eqs.  (37)  and  (38)  for  the  optimal  values 
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'  of  the 'interaction  length  £end  ar>d  the  detuning  8  for  5=0.1 
(a)  and  5=0.2  (b).  The  fact  that  decelerated  particles  move 
up  in  this  figure  corresponds  to  the  representation  of  the  elec¬ 
tron  energy  as  y=  To- Top  our  variable  y  is  proportional  to 
7(1).  A  thin  solid  horizontal  line  y-0  is  given  for  showing 
the  initial  positions  of  electrons  whose  phases  at  the  entrance 
are  uniformly  distributed  from  0  to  2tt  and  perturbation  in 
energy  is  equal  to  zero.  The  dashed  horizontal  line  in  this 
figure  corresponds  to  y=-8,  which  is  the  case  when  the 
changes  in  electron  energy  compensate  for  the  initial  mis¬ 
match  of  synchronism  between  the  wave  and  electrons.  As 
seen  in  Fig.  4,  electrons,  which  start  interaction  in  the  decel¬ 
erating  phase,  show  some  signs  of  over  saturation  to  the  end 
of  interaction  where  their  energy  y  starts  to  decrease.  Maxi¬ 
mum  changes  in  electron  energy  in  case  (b)  are  larger  that 
corresponds  to  higher  efficiencies  shown  in  Fig.  3(c)  than  in 
Fig.  3(b).  The  exponential  term  in  (37)  causes  the  speeding 
up  the  deceleration  process  in  the  final  stage  of  interaction. 
In  a  certain  sense,  this  exponential  term  leading  to  the  in¬ 
crease  of  electron  coupling  to  the  rf  field  in  the  process  of 
interaction  can  be  treated  similarly  to  the  known  tapering  of 
the  depth  of  corrugations  in  a  slow-wave  structure  resulting 
in  the  efficiency  enhancement  in  backward-wave 
oscillators.23'24  (Similar  means  of  the  efficiency  enhancement 
are  also  known  in  the  theory  of  free-electron  lasers  with  vari¬ 
able  parameters.) 

The  results  of  the  calculations  for  the  clinotron  with  a 
thick  electron  beam  are  shown  in  Fig.  5,  which  shows  the 
degradation  of  the  efficiency  rj  with  the  beam  thickness  for 
several  values  of  the  normalized  tilt  angle  a' .  The  calcula¬ 
tions  were  done  for  the  values  of  the  mismatch  of  synchro¬ 
nism  and  the  normalized  interaction  length  optimal  for  a  thin 
beam.  Recall  that  normalization  of  all  the  variables  and  pa¬ 
rameters  in  (37)  and  (39),  and  (33)— (36)  is  different  and  the 
results  shown  in  Fig.  5  correspond  to  normalization  adopted 
in  (33)-(36)  while  those  shown  in  Fig.  3  corresponded  to 
normalization  adopted  in  (37)-(39).  As  one  can  see,  when 
the  beam  thickness  reaches  its  maximum  value,  Anlax=2x0 
(here  both  A  and  x0  are  normalized  to  the  absolute  value  of 
the  transverse  wave  number  k),  the  efficiency  is  in  the  range 
of  0.5-0.6  of  its  maximum  value  for  an  ideally  thin  beam  for 
all  three  values  of  a'. 

V.  DISCUSSION 

It  should  be  instructive  to  apply  the  simple  formalism 
described  above  to  a  design  where  1 D  and  2.5D  simulations 
were  done  and  compare  the  results.  Consider  for  this  purpose 
a  recent  design  of  a  W-band  (95  GHz)  sheet-beam  traveling- 
wave  tube,10  where  ID  simulations  predicted  24%  efficiency, 
while  2.5D  PIC  simulations  yielded  the  efficiency  close  to 
22%.  That  design  was  done  for  a  planar  ridged  waveguide 
geometry  and  a  120-kV  electron  beam  of  an  elliptical  cross 
section  of  1  cm  X  0.5  mm.  The  latter  can  be  approximated  by 
a  rectangular  cross  section  of  the  same  size,  and  the  variation 
of  the  rf  field  in  the  wide  transverse  direction  will  be  ne¬ 
glected  below. 

From  the  condition  of  Cherenkov  synchronism  between 
electrons  and  the  phase  velocity  of  the  wave,  vel  ~  vph 
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FIG.  5.  Dependence  of  the  clinotron  efficiency  on  the  beam  thickness  for 
several  values  of  the  tilt  parameter  a'. 


=  a>lkz,  one  can  easily  find  the  axial  wave  number  of  the  zero 
space  harmonic  synchronous  with  the  traveling  wave:  kz 
—  ulcfi.  Correspondingly,  the  absolute  value  of  the  trans- 
verse  wave  number  of  this  harmonic,  k=(io/c)/ \ly^- 1,  for  a 
given  voltage  and  frequency  is  close  to  3.8  mm-1.  The  trans¬ 
verse  structure  of  the  axial  electric  field  in  such  a  planar 
ridged  waveguide  can  be  given  by  £.=cosh(/ar).  Corre¬ 
spondingly,  the  coefficient  C  determined  by  Eq.  (1 1)  is  equal 
to  C=(l/2)[l+sinh(K-A)//cA].  For  given  #c  =  3.8  mm"1  and 
the  beam  thickness  A =0.5  mm,  this  coefficient  is  close  to 
1.378  and  C2  — 1.9.  In  turn,  the  coefficient  D  determined  by 
(15)  is  equal  to  D=(l /2){3/4+sinh(/cA)//rA 
+sinh(2rcA)/8«-A}  that  for  given  k  and  A  yields  D—  1.99. 
Thus,  the  ratio  I=C2ID  in  this  design  is  close  to  0.955. 
Therefore,  the  transverse  nonuniformity  of  the  rf  field  should 
reduce  the  24%  efficiency  calculated  in  ID  simulations  to 
22.9%  that  is  a  little  higher  than  22%  efficiency  obtained  in 
PIC  simulations.10  This  small  discrepancy  can  be  explained, 
first,  by  our  neglect  of  the  rf  field  nonuniformity  in  the  wide 
horizontal  direction  and,  second,  by  the  fact  that  in  PIC 
simulations,  it  was  observed  that  in  the  region  of  efficient 
interaction  the  electrons  exhibit  strong  pulsations  in  vertical 
direction,  which  make  the  vertical  size  of  the  beam  in  this 
region  larger  than  initial  0.5  mm  thickness. 

VI.  SUMMARY 

A  simple  analytical  theory  is  developed  which  describes 
the  effect  of  the  transverse  nonuniformity  of  the  interaction 
space  on  the  efficiency  of  various  sources  of  microwave  ra¬ 
diation  driven  by  linear  electron  beams.  This  issue  seems  to 
be  especially  important  for  novel  sources  of  terahertz  radia¬ 
tion  where  the  beam  thickness  can  be  on  the  order  or  larger 
than  the  scale  of  transverse  nonuniformity  of  the  rf  field.  It  is 
shown  how  to  use  this  theory  for  evaluating  the  efficiency 
degradation  in  the  case  of  cylindrical  and  sheet  electron 
beams  in  various  microwave  circuits.  The  authors  hope  that 
this  theory,  while  being  unable  to  compete,  will  complement 
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well-ehborated  numerical  methods  for  detailed  design  of 
these  devices  and,  hence,  will  be  used  for  “zero-order”  esti¬ 
mates  of  the  efficiency  degradation  caused  by  the  transverse 
nonuniformity  of  the  rf  field.  Certainly,  after  minor  modifi¬ 
cations  it  can  also  be  used  for  evaluation  of  the  gain  degra¬ 
dation  in  amplifiers  with  transversely  nonuniform  interaction 
space,  since  the  gain  can  be  presented  as  G 
=  10  log{rjPb/ Pir},  where  Pb  is  the  beam  power  and  Pin  is  the 
input  power  from  a  driver. 
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The  orotron  is  a  device  in  which  a  sheet  electron  beam  excites  electromagnetic  (EM)  oscillations  in 
an  open  resonator  formed  by  a  concave  mirror  and  a  periodic  grating  plate.  The  first  spatial 
harmonic  of  a  field  near  the  grating  has  phase  velocity  close  to  electron  velocity;  therefore  electrons 
streaming  over  the  grating  surface  produce  Cherenkov  or  Smith-Purcell  radiation.  The  orotrons  can 
be  configured  for  producing  coherent  radiation  at  short  millimeter  and  submillimeter  wavelengths 
(up  to  the  terahertz  range).  Typically  the  amplitude  of  the  EM  field  in  the  orotron  is  rather  small. 
Therefore,  for  extracting  an  appreciable  part  of  electron  kinetic  energy,  it  is  necessary  to  realize 
synchronous  interaction  between  electrons  and  a  slow  wave  at  long  enough  distances.  This 
peculiarity  of  the  orotron  operation  makes  it  possible  to  develop  an  analytical  nonlinear  theory  of  the 
orotron.  Such  a  theory  showing  how  to  optimize  the  choice  of  the  interaction  length  and  the  ratio 
between  ohmic  and  diffractive  losses  is  developed  in  the  present  paper.  The  theoretical  treatment  is 
accompanied  with  consideration  of  practical  effects  limiting  the  orotron  efficiency. 
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I.  INTRODUCTION 

At  present,  there  is  a  strong  interest  in  developing  light¬ 
weight,  compact,  high-power  sources  of  terahertz  (THz)  ra¬ 
diation  for  numerous  applications.  Although  fast-wave  de¬ 
vices,  such  as  gyrotrons,  cyclotron  autoresonance  masers, 
and  free-electron  lasers  offer  very  high  power  level  at  short 
wavelengths,  they  do  not  satisfy  requirements  to  the  weight 
and  compactness.  Therefore,  the  decision  should  be  made  in 
favor  of  slow-wave  devices.  These  devices  can  be  configured 
as  amplifiers  and/or  oscillators.  Below  we  focus  on  oscillator 
configurations  only.  Among  slow-wave  oscillators  operating 
at  short  millimeter  and  submillimeter  wavelengths,  the  most 
known  are  backward-wave  oscillators  (BWOs),  whose  fre¬ 
quency  can  be  voltage  tunable  in  a  wide  range,  and  extended 
interaction  oscillators  (EIOs).  The  information  on  EIOs  can 
be  found  on  the  web  site  of  CPI,  Canada,  a  well-known 
manufacturer  of  these  devices.  The  power  level  of  BWOs 
developed  by  the  Russian  corporation  “Istok”  varies  from 
1-10  mW  at  frequencies  260-375  GHz  down  to  0.5-2  mW 
at  frequencies  exceeding  1  THz.1  This  power  level  can  be 
increased  in  the  configuration  of  BWOs  where  an  electron 
beam  propagates  at  a  small  angle  to  the  slow-wave  circuit. 
Such  configuration  is  known  as  the  clinotron;  the  latest  re¬ 
sults  of  the  development  of  clinotrons  can  be  found  in  Ref.  2. 

An  alternative  version  of  a  compact  oscillator,  which  can 
deliver  significant  power  at  wavelengths  shorter  than  1  mm, 
is  the  orotron  (this  name  is  an  acronym  for  Russian  words 
meaning  “open  resonator  with  reflecting  grating”).3  A  sche¬ 
matic  of  this  device  is  shown  in  Fig.  1 .  In  the  orotron,  elec¬ 
trons  excite  the  electromagnetic  (EM)  field  in  an  open  reso¬ 
nator  formed  by  a  smooth  concave  mirror  and  a  flat  mirror 
with  a  periodic  grating.  Due  to  grating  periodicity,  the  EM 
field  contains  spatial  harmonics.  So,  electrons  can  propagate 
in  Cherenkov  synchronism  with  the  phase  velocity  of  the 


first  harmonic  and  produce  coherent  Cherenkov  or  Smith- 
Purcell  radiation.  Typically,  the  EM  field  is  excited  in  such  a 
manner  that  the  field  of  the  zero  space  harmonics  bounces 
back  and  forth  between  the  flat  and  concave  mirror  (i.e.,  it 
has  a  small  axial  wave  number),  while  the  axial  wave  num¬ 
ber  of  the  first  spatial  harmonic  is  equal  to  2irld,  where  d  is 
the  grating  period.  Thus,  the  condition  of  Cherenkov  syn¬ 
chronism  can  be  given  as 


vet  —  vph  =  talk,  -  cd/\ .  (1) 

The  advantage  of  the  orotron  is  the  presence  of  a  high-Q 
resonator,  which  may  allow  one  to  excite  EM  oscillations 
even  at  a  relatively  low  beam  power.  Just  this  fact  had  stimu¬ 
lated  our  interest  in  this  device.  At  present,  there  are  a 
large  number  of  papers  on  the  orotron;  these  papers  were 
reviewed  in  Refs.  4-6.  Recent  experiments  with 
submillimeter-wave  orotrons7  resulted  in  delivering  up  to 
60  mW  power  at  the  370-GHz  frequency.  Theoretical 
papers  include  the  analysis  of  linear  and  nonlinear  processes 
in  one-dimensional  (ID)  and  two-dimensional  (2D)  ap¬ 
proximations  (see,  e.g.,  Ref.  6).  Of  course,  with  the  use 
of  presently  available  codes  three-dimensional  (3D)  sim¬ 
ulations  also  can  be  performed.  The  goal  of  the  present 
paper,  however,  is  completely  different.  Below  we  ex¬ 
plain  how  one  can  estimate  the  orotron  efficiency  by  very 
simple  means.  This  method  is  based  on  the  analytical  theory 
originally  developed  for  the  gyrotron.8'9  The  paper  is 
organized  as  follows.  Section  II  contains  the  formalism 
describing  the  problem,  results  of  analysis  for  an  ideal  elec¬ 
tron  beam,  and  discussion.  In  Sec.  Ill  the  effect  of  electron 
velocity  spread  is  considered.  Section  IV  summarizes  our 
paper. 
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FIG.  1.  Schematic  of  the  classical  orotron  configuration. 


harmonics  of  an  EM  filed,  can  be  described  by  two  first- 
order  equations  for  normalized  electron  energy  y=s/mc2 
and  slowly  variable  phase  d=a)t-kzz, 

$  =  -R etai/feV*},  (3) 

dz 


d-&  y 
dz!  Vy2- 1 


(4) 


II.  GENERAL  FORMALISM  AND  RESULTS  FOR  IDEAL 
ELECTRON  BEAM 

A.  Formalism 


There  are  two  facts  important  for  the  theory  developed 
below.  First,  at  submillimeter  wavelengths,  a  significant  part 
of  radiated  EM  energy  can  be  lost  in  the  resonator  walls 
having  finite  conductivity  (ohmic  losses).  Therefore,  optimi¬ 
zation  of  the  orotron  parameters  should  be  done  for  maxi¬ 
mizing  the  output  efficiency 


1 


Vd= — r  vCh 

l  +  Q 


(2) 


instead  .of  the  electronic  one,  t)el.  (In  this  regard,  optimiza¬ 
tion  of  the  orotron  has  something  in  common  with  optimiza¬ 
tion  of  gyrotrons  at  submillimeter  wavelengths.8-10)  In  (2), 
Q=Qd!  2ohm  is  the  ratio  of  the  diffractive  quality  factor  QD, 
which  accounts  for  losses  in  the  coupling  hole  shown  in  Fig. 
1 ,  to  the  ohmic  quality  factor  <2ohm;  the  total  quality  factor  in 
(2)  is  determined  by  both  kinds  of  losses,  l/Q=UQD 
+  l/Qohm-  The  diffractive  Q  factor  associated  with  the  cou¬ 
pling  hole  can  be  evaluated  as  <2D=toW/Prad,  where  W  is  the 
EM  energy  stored  in  the  resonator  and  Prad  is  the  EM  power 
propagating  via  this  hole.  When  diffractive  losses  due  to  fi¬ 
nite  sizes  of  both  mirrors  become  important,  Q0  should  ac¬ 
count  for  these  losses  as  well  (see,  e.g.,  Ref.  11).  These 
losses  are  important  for  the  selection  of  modes  with  large 
numbers  of  field  variations  on  mirrors. 

Second,  typically,  at  such  short  wavelengths,  it  is  rather 
difficult  even  to  fulfill  the  self-excitation  conditions;  thus,  the 
device  operates  at  currents  which  exceed  the  start  current 
only  slightly.  This  means  that  the  amplitude  of  EM  oscilla¬ 
tions  in  such  a  device  is  rather  small — much  smaller  than  the 
optimal  for  the  interaction  efficiency.  To  compensate  for  this 
smallness,  it  is  expedient  to  increase  the  interaction  length. 
Then,  at  least,  efficient  bunching  in  the  process  of  interaction 
may  take  place  that  should  result  in  the  efficiency  enhance¬ 
ment  in  comparison  with  a  device  of  a  shorter  length. 

Let  us  illustrate  the  consequences  from  this  fact  (small 
EM  amplitude,  long  interaction  length)  with  the  use  of  equa¬ 
tions  for  electron  motion  and  the  balance  equation  for  the 
EM  power  produced  by  an  electron  beam  and  the  power  of 
EM  losses  in  a  cavity.  For  simplicity,  we  will  ignore  the 
electron  velocity  spread  and  the  transverse  nonuniformity  of 
the  interaction  space  (these  factors  will  be  discussed  later). 
As  shown  elsewhere,  '  ID  electron  motion  in  a  resonator, 
where  electrons  are  in  synchronism  with  one  of  the  spatial 


In  (3)  and  (4),  we  introduced  a  number  of  normalized  param¬ 
eters:  the  normalized  axial  coordinate  z'  =  ojz/c,  the  normal¬ 
ized  amplitude  of  the  first  spatial  harmonic  of  the  EM  field 
ai=eAl/mco)  (the  function  f(z')  describes  its  axial  struc¬ 
ture),  and  the  normalized  axial  wave  number  h=kzcl to, 
which  is  the  inverse  phase  velocity  of  the  first  spatial  har¬ 
monic  normalized  to  the  speed  of  light.  Also,  the  first  term  in 
the  right-hand  side  (RHS)  of  (4),  1 1  (i  (fi=vellc  is  the  elec¬ 
tron  velocity  normalized  to  the  speed  of  light),  was  written 
with  the  account  for  the  general  relation  between  the  electron 
velocity  and  energy,  y=  1  /  vl  -  01.  The  boundary  conditions 
at  the  entrance  are  y{z'=  0)  =  y0,  r)(z'=0)  =  i‘)0,  where  initial 
phase  #0  is  uniformly  distributed  from  0  to  2m 

In  these  notations,  the  electronic  efficiency  can  be  deter¬ 
mined  as 


Here  L'  =  u>Llc  is  the  normalized  interaction  length.  The  bal¬ 
ance  equation  can  be  given  as 

£ W=  Vcpb-  (6) 

Here  Pb=h^b  is  the  beam  power.  Beam  voltage,  Vb,  and 
the  normalized  initial  energy  of  electrons,  y0,  are  related  as 
y0=  1  +eVb/mc2.  Since  the  amplitude  of  the  zero  spatial  har¬ 
monic,  A0,  is  typically  much  larger  than  amplitudes  of  other 
harmonics,  the  EM  energy  stored  in  the  resonator  volume  V 
can  be  determined  as  W—{a\I%tt)V. 

When  perturbations  caused  by  the  EM  field  in  both  the 
electron  energy  and  phase  are  small,  one  can  readily  get  from 
(3)-(6)  the  starting  conditions,  which  determine  the  start  cur¬ 
rent  of  the  device  (such  current  was  analyzed  elsewhere13). 
Below  we  consider  the  situation  when  the  beam  current 
slightly  exceeds  its  starting  value.  In  such  a  case  the  oscilla¬ 
tions  are  excited  with  small  amplitude.  Hence,  in  (4),  we  can 
treat  changes  in  the  electron  energy  in  the  process  of  inter¬ 
action  as  perturbations  (y=  y0-  y,  where  |y|  y0),  while  the 

changes  in  phase  can  be  significant  at  a  long  enough  dis¬ 
tance.  Under  these  assumptions,  Eqs.  (3)  and  (4)  can  be  re¬ 
written  as 

~  =  Re(or1/(z')e'fl),  (7a) 

dz 
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FIG.  2.  Electron  bunching  efficiency  t]h  and  optimal  detuning  5  as  functions 
of  the  normalized  interaction  length  <;0UI. 


M  1  y 

^=J0 


(7b) 


and,  hence,  these  two  equations  can  be  reduced  to  one  non¬ 
linear  pendulum-like  equation, 

•p,  Q.  1 

:Re{«^(zV*}.  (8) 


dz 


t 2  ' 


(To-  1) 


3/2 


The  amplitude  of  the  first  spatial  harmonic  o,  as  well  as  the 
function  f(z')  describing  its  axial  structure  can  be  assumed 
real.  Then,  introducing  a  new  normalized  axial  coordinate 
we  can  reduce  (8)  to 

rp,  Q. 

77  =/(<?)  Re(/8).  (9) 

<?<; 


This  equation  should  be  supplemented  with  the  boundary 
conditions  at  s=0  for  the  electron  phase  (i7(0)  =  -il0 
e[0;27r])  and  its  derivative  <?i7/<?q|,=0=<5,  where  <5 
=  (l/y3o-/i)/V«1/(^-l)372  is  the  normalized  initial  detun¬ 
ing,  which  is  the  same  for  all  electrons  in  an  ideal  (zero 
velocity  spread)  electron  beam. 

Equation  (7b)  shows  that  the  changes  in  electron  energy 
are  proportional  to  the  phase  derivative.  Therefore,  the  elec¬ 
tronic  efficiency  (5)  can  be  rewritten  as 


Ver 


(ro+1)3 

L  7o-l  J 


1/4  r  1 

van 


j_r2w 

2ttJ0 


dd_ 

<?S 


(10) 


The  expression  in  figure  brackets  characterizes  the  electron 
bunching.  Therefore  it  will  be  called  below  the  bunching 
efficiency  and  denoted  by  r/b.  As  shown  elsewhere,914  its 
value  can  be  larger  than  1.  This  efficiency  optimized  with 
respect  to  the  normalized  detuning  is  shown  as  the  function 
of  the  normalized  interaction  length  for  the  device  with  a 
Gaussian  axial  structure  of  the  EM  field  in  Fig.  2  reproduced 
from  Ref.  14.  (Here,  the  Gaussian  distribution  corresponds  to 
the  case  when  at  the  entrance  and  the  exit  the  field  amplitude 
is  e  times  smaller  than  in  the  middle  of  the  resonator.)  It 
should  be  emphasized  that  because  of  increasing  the  interac¬ 
tion  length,  the  electronic  efficiency  (10)  scales  at  small  am¬ 
plitudes  as  \fa[,  while  in  devices  of  a  fixed  length  it  scales  as 
(see,  e.g.,  Ref.  12).  The  expression  in  square  brackets  in 
(10),  at  low  voltages,  reduces  to  [(y(j+1)3/(y0-l)]1/4 


~(8 mc2leVb)w.  Note  that,  as  the  voltage  reduces,  the  EM 
field  amplitude  reduces  as  well,  as  follows  from  the  balance 
equation. 

In  accordance  with  the  steps  made,  the  balance  equation 
(6)  can  be  rewritten  as 


1  +  2 


Vb_ 

‘'out 


(11) 


Here  the  normalized  current  parameter  /  is  equal  to 


2  ej\2aZwv 
me  a0  H 


In  (12),  j  is  the  electron  current  density  in  the  interaction 
space,  a0  is  the  amplitude  of  the  zero  spatial  harmonic  of  the 
EM  field  normalized  in  the  same  fashion  as  ax  above.  The 
last  ratio  wv/H  is  the  ratio  of  the  thickness  of  a  sheet  electron 
beam  in  the  vertical  direction,  w„,  to  the  distance  H  between 
two  mirrors  shown  in  Fig.  1.  In  transforming  (6)  into  (11)  we 
took  into  account  not  only  normalization  of  a  new  variable 
gout,  but  also  the  fact  that  the  resonator  volume  increases  with 
the  length,  as  the  amplitude  of  EM  oscillations  gets  smaller. 
It  was  also  assumed  that  the  beam  horizontal  width  is  equal 
to  the  resonator  size  in  the  direction  perpendicular  to  the 
schematic  shown  in  Fig.  1 . 

So,  Eqs.  (2)  and  (9)— (1 1)  form  a  set  of  equations,  which 
can  be  used  for  optimizing  the  Q  ratio  Q  and  the  normalized 
interaction  length  %ut  for  maximum  output  efficiency. 


B.  Results  and  discussion 

Optimization  with  respect  to  the  Q  ratio  readily  yields 
the  optimum  value 

2oP,=  l/2.  (13) 

This  value,  obtained  with  the  account  for  the  dependence  of 
the  EM  field  amplitude  a,  on  the  Q  ratio  given  by  (11), 
shows  that  1/3  of  the  radiated  microwave  power  should  be 
lost  in  the  resonator  walls,  while  2/3  can  be  extracted  via  a 
coupling  hole  shown  in  Fig.  1. 

Optimization  of  the  normalized  length  should  be  done 
with  the  use  of  Eqs.  (2),  (10),  and  (11).  From  (1 1)  it  follows 
that 


da , 


2 


^%ut  1  +  Q  ^  out  '  ^out 


(14) 


while  optimization  of  the  output  efficiency  with  respect  to 
this  length  yields 


1—dVb  1  dai  r, 
V«i- — +  TT=Vb- —  =  0. 


(15) 


*out  2  yo^  ^  out 

With  the  use  of  (14)  and  (11),  Eq.  (15)  can  be  reduced 


to 
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*  (— j*)  =\-  (16) 

\  Vb  flyout '  opt  1 

For  the  case  of  the  Gaussian  axial  structure  of  the  EM  field, 
solution  of  (16)  can  be  found  with  the  use  of  Fig.  2.  It  is 
%ut, opt- 2-65;  the  corresponding  value  of  the  bunching  effi¬ 
ciency  is  equal  to  1.78.  As  is  seen  in  Fig.  2,  this  normalized 
length  is  smaller  than  the  one  optimal  for  the  bunching  effi¬ 
ciency. 

For  these  optimal  values  of  the  Q  ratio  and  normalized 
interaction  length,  the  normalized  amplitude  of  the  first  spa¬ 
tial  harmonic  of  the  EM  field  is  equal  to  al  opt  —  0.2247,  and 
the  maximum  output  efficiency  is  equal  to 

%u.,max-0.567'/2.  (17) 

So,  the  results  of  our  optimization  show  that  the  maxi¬ 
mum  output  efficiency  depends  on  the  normalized  beam  cur¬ 
rent  parameter  defined  by  (12).  Most  of  the  parameters  in 
(12)  are  essentially  the  same  as  those  analyzed  in  Refs.  13 
and  15,  which  allows  us  to  make  use  of  the  analysis  done 
there.  In  particular,  it  was  found15  that  the  maximum  ratio  of 
the  amplitude  of  the  first  spatial  harmonic  to  that  of  the  ze¬ 
roth  one  can  be  achieved  when  the  grooves  shown  in  Fig.  1 
are  a  quarter-wavelength  deep  and  their  width  is  equal  to  half 
of  the  period  d.  In  such  a  case,  the  ratio  a{/a0  in  (12)  is 
equal  to  2/ir.  Also,  the  ohmic  Q  factor  in  such  a  case,  as 
shown  in,15  is  equal  to 

H 

Gohm=  /  ,  \  •  (18) 


Here  Ssk  is  the  skin  depth,  which  we,  for  practical  reasons, 
will  take  equal  to  twice  its  theoretical  value 
=  (1/2it)  'iKc/a  (cr  is  the  wall  conductivity).  Equation  (18)  is 
very  important.  In  the  absence  of  grooves,  the  ohmic  Q  fac¬ 
tor  of  a  resonator  formed  by  two  smooth  plates  is  equal  to 
Qohm=HI2Ssk.  Comparison  of  this  formula  with  (18)  shows 
that  the  presence  of  grooves  reduces  the  ohmic  Q  by  the 
factor  of  (3  +  X/t/)/2.  Once  we  take  into  account  that,  in 
accordance  with  the  Cherenkov  synchronism  condition  (1), 
\ld—\lpQ,  we  come  to  the  conclusion  that  just  a  large  num¬ 
ber  of  deep  and  narrow  grooves  forming  a  grating  plate  plays 
the  crucial  role  in  the  enhancement  of  ohmic  losses.  So,  be¬ 
low  we  will  determine  the  ohmic  Q  factor  as 


_  H  A, 

4*  1  +  3A) 


(18a) 


because  (18a)  better  illustrates  the  consequences  of  lowering 
the  beam  voltage  for  the  device  performance.  At  low 
voltages,  the  normalized  electron  velocity  is  equal  to 
'j2eVblmc2—  \jVb(kV)1 16. 

Lastly,  the  thickness  of  a  beam  in  the  vertical  direction, 
Wy,  should  be  small  enough  for  avoiding  efficiency  degrada¬ 
tion  because  of  the  vertical  nonuniformity  of  the  slow  wave 
localized  near  the  grating.  As  shown  in  Ref.  12,  this  degra¬ 
dation  becomes  important  (the  efficiency  degrades  by  a  fac¬ 
tor  of  2)  when  the  beam  thickness  is  on  the  order  of  \/30/7 r. 
So,  expressing  the  distance  between  mirrors,  77,  in  terms  of 


half  wavelengths,  H—q(\/2),  where  index  q  shows  the 
number  of  field  variations  in  the  vertical  direction,  the  ratio 
wvIH  in  (12)  can  be  estimated  as  2/30/ ttcj.  With  all  this  in 
mind,  the  expression  for  the  normalized  current  parameter 
(12)  can  be  rewritten  as 


7  = 


8  X  ejfk 


Pi 


7 r  8sk  me3  1  +  3(30 


(19) 


This  equation  shows  that  to  maintain  the  value  of  this 
parameter  fixed  as  the  wavelength  shortens  the  beam  current 
density  should  increase  inversely  proportional  to  X5/2.  Also, 
it  follows  from  (19)  that  at  low  voltages  (when  3/30<§  1)  the 
beam  current  parameter  7  is  linearly  proportional  to  the  beam 
voltage. 

To  illustrate  the  meaning  of  this  equation,  let  us  consider 
an  example:  the  wavelength  equal  to  0.8  mm  (at  this  wave¬ 
length,  the  skin  depth  for  copper  is  close  to  0.22- 10~4  cm), 
the  beam  voltage  equal  to  4  kV,  and  the  beam  current  den¬ 
sity  equal  to  30  A/cm2.  For  these  numbers,  the  normalized 
beam  current  given  by  (19)  is  equal  to  0.38  - 10-4.  (The  thick¬ 
ness  of  such  a  beam  in  vertical  direction  should  not  exceed 
32  yttm.)  This  value  of  the  normalized  beam  current  corre¬ 
sponds  to  the  normalized  amplitude  of  the  first  spatial  har¬ 
monic  of  the  EM  field  equal  to  0.85- 10“5.  For  a  given  volt¬ 
age  and  this  amplitude,  the  optimal  normalized  length 
s0m,opt=2.65  corresponds  to  the  ratio  of  the  interaction  length 
to  the  wavelength  L/X  — 6.4.  For  such  small  amplitude,  the 
maximum  output  efficiency  given  by  (17)  is  close  to  0.35% 
only.  Of  course,  the  total  radiated  microwave  power  in¬ 
creases  with  the  width  of  a  system  in  a  wide  horizontal  di¬ 
rection.  For  instance,  when  this  width  is  equal  to  3  mm,  the 
total  beam  current  is  about  30  mA.  Correspondingly,  the 
beam  power  is  120  W  and  the  power  of  outgoing  radiation  is 
0.42  W.  Recall  that  here  we  consider  the  beam  current  den¬ 
sity  in  the  interaction  space,  which  can  be  increased  by  com¬ 
pressing  a  beam  in  the  transition  region  between  electron  gun 
and  the  interaction  space,  when  such  a  region  exists.  For 
instance,  if  the  current  density  is  increased  due  to  beam  com¬ 
pression  from  30  A /cm2  taken  above  to  300  A /cm2,  the 
beam  current  parameter  7  (19)  also  becomes  10  times  larger, 
which  increases  the  output  efficiency  (17)  from  0.35%  to 
about  1.1%.  Correspondingly,  the  total  power  of  radiation 
increases  from  0.42  to  13  W.  (Now  we  assume  that  the  total 
beam  current  is  0.3  A.) 

The  beam  compression  can,  in  particular,  be  realized  in 
miniaturized  plasma-cathode  electron  guns16  suitable  for 
short-wavelength  applications.  Originally,  similar  guns17 
were  developed  for  plasma-assisted  microwave  oscillators.18 
In  such  guns,  the  beam  compression  can  be  realized  by 
proper  axial  profiling  of  the  neutral  gas  density,  because  im¬ 
pact  ionization  of  these  neutrals  by  a  propagating  electron 
beam  produces  the  ions  providing  ion  focusing  of  the  beam. 


III.  EFFECT  OF  ELECTRON  VELOCITY  SPREAD 

The  method  proposed  has  some  limitations  once  the 
spread  in  electron  velocities  is  taken  into  account.  Let  us 
come  back  to  the  equation  for  the  slowly  variable  phase  (4) 
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*  arid  determine  the  limitation  on  the  interaction  length  caused 
by  electron  velocity  spread. 

The  equation  for  a  slowly  variable  phase  (4)  can  be  re¬ 
written  as 


JL_ 

dZ’~  PpH 


(20) 


where  /3ph  is  the  phase  velocity  of  the  first  spatial  harmonic 
of  the  wave  normalized  to  the  speed  of  light.  Above  we  were 
focused  on  the  effect  of  changes  in  the  electron  velocity  and 
energy  in  the  process  of  interaction  on  the  phase  bunching. 
Now,  let  us  evaluate  the  effect  of  the  spread  in  initial  veloci¬ 
ties  on  the  synchronism  between  electrons  and  the  wave.  In 
accordance  with  (20),  electrons  having  initial  normalized  ve¬ 
locity  different  from  the  mean  value  /30,  say  /j0=/30+ A/30, 
experience  additional  phase  shift,  which  can  be  evaluated  as 
A'd=(-l/yS0)(A^0/yS0)z'.  Assume  that  the  velocity  spread  is 
acceptable  when  all  electrons  stay  in  approximately  the  same 
phase  with  respect  to  the  wave,  i.e.,  when  the  phase  shift 
caused  by  the  velocity  spread  does  not  exceed  it/2.  This 
limitation  immediately  gives  us  the  maximum  ratio  of  the 
interaction  length  to  the  wavelength 


^max  A)  (')]} 

X  =4  sp0’ 

where  S/30= Aj80//30  is  the  relative  spread,  which  can  be 
given  in  percents.  Taking  into  account  the  Cherenkov  syn¬ 
chronism  condition  (1),  Eq.  (21)  can  be  rewritten  in  terms  of 
the  number  of  grating  periods  N=L/d  as 

Nma=\/48(30.  (21a) 


obeys  restrictions  (21)  and  (21a)  and  whose  grating  period  is 
constant,  the  electron  bunch  is  formed.  Then,  deceleration  of 
this  bunch  takes  place  in  the  second  stage  where  the  grating 
period  varies  with  distance.  This  concept  of  a  device  with 
tapered  parameters  is  known  in  traveling-wave  tubes19  and  in 
free-electron  lasers.20,21  The  wall-plug  efficiency  of  such 
low-efficiency  devices,  of  course,  can  also  be  greatly  in¬ 
creased  by  using  depressed  collectors.  This  means,  however, 
complicates  the  device  and  increases  its  weight  and  sizes. 


IV.  SUMMARY 

A  simple  theory  presented  in  this  paper  allows  develop¬ 
ers  to  easily  evaluate  the  efficiency  and  the  most  important 
constraints  in  the  design  of  the  orotrons  intended  for  opera¬ 
tion  at  submillimeter  wavelengths,  up  to  the  THz  range.  As 
follows  from  the  results  presented,  two  of  the  most  critical 
issues  for  operation  at  such  high  frequencies  are  the  ohmic 
losses  in  small-period  gratings  and  relatively  low  electron 
current  density.  To  solve  the  first  problem,  a  variety  of  alter¬ 
native  concepts  of  periodic  structures  can  be  proposed  and 
analyzed.  Manufacturing  of  such  structures,  when  necessary, 
can  be  done  with  the  use  of  novel  methods  of  microfabrica¬ 
tion  (see,  e.g.,  Ref.  22).  Thus,  the  structures,  which  seemed 
unrealistic  some  time  ago,  can  be  revisited  today.  To  solve 
the  problem  of  low  current  density,  new  emitters  with  high 
current  density  as  well  as  methods  of  beam  compression  re¬ 
quiring  a  careful  design  of  electron  optical  systems  with  the 
use  of  presently  available  numerical  codes,  should  be  used. 
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25.  Coming  back  to  the  optimal  ratio  L/X  —  6A  found  in  the 
example  considered  above,  we  can  conclude  that  for  realiz¬ 
ing  this  ratio  the  velocity  spread  should  not  exceed  0.5%. 

This  issue  can  be  analyzed  in  a  more  general  manner 
with  the  use  of  the  formalism  developed  above.  Combining 
the  optimal  value  of  the  normalized  amplitude  of  the  first 
harmonic  of  the  EM  field,  a,  opt  —  0.224/,  with  the  optimal 
value  of  the  normalized  length  sout,opt=2-65  and  the  limita¬ 
tion  on  the  interaction  length  given  by  (21)  or  (21a),  one  can 
readily  find  that  to  fulfill  this  limitation  the  electron  current 
density  in  the  interaction  space  should  obey  the  following 
condition: 


k  ej X2  A, 

8 A  me3  !  +  3/30 


>  155(<5/30)2. 


(22) 


When  this  condition  cannot  be  fulfilled,  the  interaction 
length  should  be  chosen  in  accordance  with  (21),  while  the 
nonnalized  amplitude  will  be  determined  by  the  beam  cur¬ 
rent,  in  accordance  with  the  balance  equation  (11).  As  a  re¬ 
sult,  the  normalized  length  sout  will  be  smaller  than  its  opti¬ 
mal  value  found  above  and,  correspondingly,  the  output 
efficiency  will  be  smaller  than  its  optimal  value. 

To  enhance  the  interaction  efficiency,  one  can  consider  a 
two-stage  device,  where  in  the  first  stage,  whose  length 
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